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and for overall project management. Dr. T. L. Wong 
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for the experimental aspects of the study. 

This report has been assigned NASA report number CR- 11 5087 
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ABSTRACT 


This report describes results from a program to investigate the 
stabilizing influence of acoustic cavities on acoustic modes of 
combustion Instability*. The primary objective of the program 
was to develop sufficient understanding, and the necessary data, 
concerning cavity behavior so that these cavities can be effectively 
designed for combustion stabilization. 

During this program, several areas of investigation were pursued: 
Increased analytical capability was developed by removing or 
improving some of the previously_used approximations and allow- 
ing for additional effects and configurations. Test firings 
were made to experimentally define an optimum slot width (slot 
width corresponding to maximum damping); modifications have been 
made to the analytical model to improve the ability to predict the 
optimum sloil width. Also, test firings have been made to investi- 
gate the effects of cavity location and of cavity multiplicity 
(multiple rows of resonators) on stability. Detailed cavity 
temperature measurements have been made. Analytical and sub- scale 
experimental studies of unconventional cavity configurations (which 
emphasize hardware constraints) have been made; full-scale stability 
tests with these cavities have been made. Eullescale tests were 
also made to investigate the. effect of varyin g the film coolant flow- 
rate and nominal operating conditions on cavity stabilization. 
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Results from all phases of the program provide a substantial basis 
for. .cavity design. The utility of these cavities has_again been 
demonstrated. Dynamic stability of the engine was obtained with a 
wide range of cavity configurations and conditions. 
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INTRODUCTION AND SUMMARY 


Acoustic modes of combustion instability in a rocket engine can be prevented 
through the use of acoustic absorbers. These absorbers are comprised of 
acoustic resonators distributed in some manner along the interior walls of 
a thrust chamber. One such arrangement, which is particularly attractive 
from a_design and manufacturing standpoint, is a single row^o£_acoustic 
resonators along the periphery of the injector. The term "acoustic cavity" 
has been loosely used to describe this simple arrangement of resonators, 
generally quarter wave resonators. A typical acoustic cavity arrangement 
is shown schematically in Fig. 1. These cavities are simply narrow slots, 

' either axially or radially directed, with uniform cross-sectional area. The 
slots are partitioned to prevent circumferential flow of hot gases. How- 
ever, interest is not restricted to these simple slots. Any fom of acoustic 
resonator should be useful, although seme forms may be more effective than 
others. 

This report describes results from a program to investigate the stabilising 
influence of acoustic cavities. The primary objective of the program was to 
develop sufficient understanding, and the necessary data, concerning cavity 
behavior so that these cavities can be effectively designed for combustion 
stabilization. 

The current investigation is based heavily on results from, and is- largely 
a continuation of a previous study performed as part of the 121 Ascent Engine 
program, which is described in Ref. 1 and 2. During that study substantial- 
progress was made concerning the influence of cavities on stability and 
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Typical Cavity Configuration 


how to desig n them. Furt her, the stabilization effectiveness of acoustic 
cavities was demonstrated. An analytical model was developed to predict 
the damping provided by a cavity. Also, subscale tests were made with a 
combustion-driven oscillator, a T-bumer, to measure the high-amplitude 
acoustic characteristics of the cavity. In addition, a number of full- 
scale hot firings were made with LM-ascent-engine-type hardware and with an 
unbaffled injector. During this previous program, dynamic stability was 
obtained with five different cavity arrangements while instability was 
easily triggered without the cavity. 

During the current program, several areas of investigation have been pursued: 
Increased analytical capability was developed by removing or improving some 
of the previous!y_used approximations and allowing for additional effects 
and configurations . Test firings were made to experimentally define an 
optimum slot width; modifications have been made to the analytical model to 
improve the ability to predict the optimum slot width »_ Also, test firings 
hav e been made to investigate the effects of cavity location and of cavity 
multiplicity (multiple rows of resonators) on stability. Detailed cavity 
temperature measurements have been made. Analytical and sub-scale experi- 
mental studies of unconventional cavity configurations (which emphasize 
hardware constraints) have been made; full-scale stability tests with these 
cavities have been made. Full-scale tests were also made to Investigate 
the effect of varying the film coolant flowrate and nominal operating con- 
ditions on cavity stabilization. 


For convenience the program was divided, into three phases: Phase I - Develop 

Design Criteria, was comprised of an analytical study, a subscale experimental 
study, and a full-scale experimental study. Phase II - Unconventional Cavity 
Testing, Involved a full scale test series (originally a series of tests to 

verify design procedures was planned but this series was later deleted In 

favor of the unconventional- cavity series). Phase III - Film Cooling and 
Off-Nominal Testing, involved additional full scale testing to evaluate the 
indicated- effects. - - — — 

SUMMARY OF RESULTS 

The analytical effort was directed toward strengthening the analytical model 
for cavity damping developed during the previous program. The principal 
weakness of the model was the inability to predict the observed cavity width 
(open area) corresponding to maximum damping or stability, .iajor improve- 
ments have been made in the accuracy with which the equations comprising 
the model_are solved. Approximate criteria were developed to assure con- 
vergence of series expressions appearing in the model. Uncertainties rela- 
tive to aperture end corrections and average impedance values were removed. 

The model was extended to specifically include radially directed cavities. 

The effects of combustion driving and various loss processes were simulated to 
same extent by introducing gain/loss type boundary conditions. The most 
significant improvement in aecuracy was obtained by introducing an iterative- 
approximation. These improvements allowed the model equations to be accurately 
solved. 
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The model was found to properly predict the observed cavity width corresponding 
to maximum damping or stability If a cavity- Impedance amplitude parameter 
rp/>p 0 was chosen large enough. This corresponds to the amplitude for which 
the cavity is designed and is considered to be a satisfactory explanation of 
observed behavior. As currently developed, this model is considered well suited 
for the design of all kinds of acoustic absorbers. It forms to basis of the 
recommended design approach. 

Moreover, analytical studies were made of unconventional cavity configurations. 

An analytical model was developed for a rather general cavity configuration. 

It was demonstrated that an "L" shaped resonato. exhibits bothHelmholtzand 

quarterwave-type resonances which allows a single cavity to be tuned within 
limits to damp two modes of instability. Also , it was shown that a quarter- 
wave resonator occupies a minimum volume relative to Helmholtz and intermediate 
resonators with the same open area and resonant frequency. In many cases, 
unconventional cavities offer significant advantages over conventional configura- 
tions. 

Subscale tests were made, using the T-bumer to investigate cavity damping 
and impedance with (l) multiply tuned cavities and (2) unconventional cavities. 

The results were not entirely satisfactory but do provide some support for the 
analytical formulation. The data indicate that multiply tuned cavities can be 
used to increase damping bandwidth. In addition, needed cavity impedance information 
was obtained from the unconventional cavity tests. 
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Pour different series of full-scale test firing were made to investigate various 
cavity effects on stability. During the Phase I full scale testing, the in- 
fluence of cavity width, axial location and multiplicity on engine stability 
were measured. The results indicate a relatively large range of cavity dimen- 
sions exists which will effect stability. No slot width up to 1.0 Inches was 
found "too" wide to stabilize the engine; the maximum effective slot width 
was not reached during the program. In addition, the engine was readily 
stabilized with radially directed slots located up to 1.0 inches (midwidth) 
downstream from the injector face; thus, _the stabilization ability of the slot 
is not highly sensitive to axial position as seme previous data had suggested. 
Moreover, the engine was also stabilized with two and three slot combinations; 
however, the only evident advantage of multiple slots over a simple slot with 
equivalent open area is to introduce multiple resonant-frequencies of the 
cavities. Furthermore, during these tests the influence of engine operating 
conditions was removed by precise control of the engine mixture ratio and 
chamber pressure. 

The influence of operating conditions on cavity stabilization was investigated 
during the off nominal test series. The test results exhibit a gradual and 
consistent variation of stability with mixture ratio and chamber pressure. 

These results probably reflect simply a variation in the instability driving 
processes rather than a cavity effect. The stabilization ability of a cavity 
does not appear highly sensitive to engine operating conditions. 
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The influence of the film coolant flow rate on cavity stabilization was also 
investigated. The stability results indicate little influence, if any, of 
film coolant flowrate over the range investigated. The observed weak dependence 
of cavity stabilization on film coolant flowrate removes concern for a degrada- 
tion of stability accompanying increased film coolant flowrate which has been 
observed with acoustic, liners. 

The stabilization ability of four different unconventional. cavity configura- 
tions (two "L" shaped and two intermediate configurations) was -evaluated. 

All .four produced dynamic stability. This demonstrated effectiveness allows 
greater design flexibility to be introduced and more efficient utilization, of 
available space. Moreover, some stability gains were achieved} the 0.3-inch- 
wide intermediate cavity stabilized the engine whereas none of the 0.3- inch- 
wide conventional cavities completely stabilized the engine. Furthermore, 
the unconventional cavities were designed for two widely separated cavity 
temperatures and, thereby, resonant frequencies. Nevertheless, all of them 
stabilized the engine, a result which suggests a substantial stability margin. 

Relatively detailed cavity temperature measurements were made during essentially 
all hot firings. The resultant temperature data appear more precise and con- 
sistent than any previously available. These data can be averaged with 
greater confidence than previously. Moreover, the influence on cavity tempera- 
ture of cavity width and depth, and of engine operating conditions and film 
coolant flowrate are reasonably clear from the data. Further, the data-frem 
the unconventional cavities su gg est the temperatures therein are largely 
equivalent to those in conventional cavities. 
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During this program a low-cost bomb was developed which exhibits the same 
yield as previous bombs.. but with, a substantial cost saving. 

Results from all phases of the program provide a substantial basis for cavity 
design. The utility of. these_cavities has again been demonstrated. 
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ANALYTICAL STUDIES 


An analytical approach, described in Ref. 1, for the analysis and design, 
of acoustic absorbers was developed .during the LM ascent engine acoustic 
cavity study. This approach involves analytically estimating a damping 
contribution provided by the absorber, the contribution being expressed as a 
temporal damping coefficient. During the current program, analytical studies 
were performed to strengthen the previously developed approach and to analyze 
contributory processes relating to cavity damping. Hie approach is based 
on the use of an approximate separation of contributory processes and, 
further, of a quasilinear representation of the oscillatory behavior so that 
the analytical methods of acoustic theory may be employed . With these 
approximations a relativel y tractab le analytical representation of cavity 
effects is obtained. Although it may be desirable to eliminate these 
approximations in favor of a more nearly exact approach, a more nearly 
complete nonlinear analytical model is like ly to be unwieldy a nd intractable; 
therefore, its utility would be restricted. 

Calculation of a damping coefficient requires analytical representations 
for both the-jwave motion_in the combustion chamber and that in the cavity 
itself. 

ANALYTICAL MODEL 


The analytical model developed to describe the stabilizing effect of the 
cavity was based on the calculation of a temporal damping coefficient con- 
tribution arising from the absorber. This parameter is believed to most 
strongly reflect the stabilizing effect of the cavities. The damping coeffi- 
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cient was chosen because, if a system of equations describing the low- 
amplitude (linear) wave motion in a combustion chamber including the effects 
of combustion was solved, a time dependence would be obtained for each made 
of the fora: 

e~ cos 


The overall damping coefficient is « N and the oscillatory frequency is 
ujj/2* . If is negative, the system is unstable and the amplitude grows 
exponentially. Furthermore, <* N is approximately the sum of a series of 
contributions, i.e., 

°Vi V (i) 

XI v ' 


where driving processes contribute negative a^'s and damping processes 
contribute positive values. Equation 1 was developed and explained more 
fully in Ref*_lj it is also discussed by Hart (Ref. 3 ) * Ihe utility of 
this approach is that a good estimate of particular contributions to this 
sum, e.g., a contribution due to an acoustic cavity, is obtained from an 
analysis which neglects other gains and losses. Consequently the analytical 
model was based on the calculation of a damping coefficient contribution, 
an due to the cavity, while ignoring other gain and loss processes. 
These other processes are, generally, too poorly understood to justify 
attempts to calculate the remaining a^’s. Therefore, it is assumed that., 
the overall «„ will be positive and the engine stable if the contribution— 
due to the cavities is made "large enough. 11 The damping coefficient due to 
the acoustic cavities is believed to be the best available measure of the: 


stabilizing effect of the cavities, even under high amplitude (nonlinear) 
conditions. 
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The damping coefficient may be suitably estimated by a calculation that 
neglects other gain and loss effects, an approach which is widely used In 
studies of solid-propellant instability. Toward this end, it is appropriate 
to solve the wave equation for a chamber of the size and shape of the com- 
bustion chamber. This solution is complicated somewhat by the fact that 
the boundary condition, the wall impedance is nonuniform; therefore, the 
equation cannot be solved by the usual separation of variables. A suitable 
alternate technique is to convert the wave equation and boundary conditions 
to an integral equation and obtain an approximate solution to the latter 
equation by a variational techni que, as suggested by Morse and Ingard 
(Ref. 4, PS* 56l). The integral equation can be solved by other techniques, 
but the variational method is probably the most suitable. 


The damping coefficient is obtained by solving the wave equation: 



J_ if? 

c 2 3t 2 


( 2 ) 


subject to the boundary conditions: 



-j ~ at cavity 
0. elsewhere 



where JS is the oscillatory pressure, c is main-chamber sound velocity* c, 
is the specific impedance of the cavity, and (3 is the complex angular 
frequency. The analysis applies equally well to ali-kimde of resonators 
if an approproprinte impedance expression is used. 
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Equations 2 and 2 may be rewritten as an integral equation, 

P(r) * ^G(r [r^ ^V o p(r o ) dS o (4) 

where the integral is a surface integral and G(r|r 0 ) is a Green's function, 
which satisfies a zero-gradient boundary condition. The Green's function 
also satisfies the differential equation 

+ k 2 G = - 6(?-? 0 ) (5) 


where k = (3/c and 6 (r-r Q ) is a Dirac delta function. Morse (Ref. 5$ pg. 791-03^) 
describes several ways of obtaining Green's functions. One suitable form is 




( 6 ) 


where 
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An alternate derivation of the integral equation ma y be obtained by employing 
a series of separated solutions. This alternate approach is less general and 
flexible than the former but has the advantage of not requiring the introduc- 
tion of a Green's function at the outset. The_approach may be illustrated 

for a cylinder with rigid walls except for one end (z ® L) which has an 
arbitrary boundary condition. The separation-of-variables technique is used 
without initially applying the boundary condition for the nonrigid end* A 
sum of the resultant solutions is then used in an effort to satisfy the re- 
maining boundary condition. Thus, a solution is sought of the form 

p = V* (a cos m0 + b sin m0)cos k z J (a r/r ) (7) 

* mn mn mn m ®n " 

m,n 

where a mn and b mn are undetermined coefficients. 


The pressure gradient at z = L may be obtained by differentiation, i.e.. 



- V (a cos m0 + b sin m0)k sin k L ) 

/ . mn mn mn mn m mn w 

m,n 


( 8 ) 


By treating this gradient as known, the coefficients 
tained from orthogonality of the functions, thus 

Jfr * I 2=l cos m6 o J m ^mn r o/ r v )de o r o dr o 



and nay bs ob- 


(9) 


The coefficient b^ is similar but with cos m 0 replaced by sin m e Q . 
Equation (7) then becomes 
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This equation is identical with the equation obtained by inserting Eq. (6) 
into Eq* (4) (with z < z Q and z 0 L) . Essentially, this approach is 
simply an alternate way of defining a Green 1 a_function. However, it may be 
more appealing to individuals who are unfamiliar with Green's functions. 


Through application of a variational procedure suggested by Horse (Ref* 4, 
Pg*5.6l)> a characteristic equation was obtained from Eq. 4 , i.e., 

y ^mnv + mnv ^ r l r <>) ^ * 0 OX)' 

where y ° l/C and an approximate pressure distribution at the slot (eigen- 
function) of the form A i — has been employed for the variation, the value 

mnv 9 

of A being optimized by the variation. 


For an axially directed slot located in the injector face, with the pressure 

distribution for the no-slot case used for t — and ajuniform cavity 

mnv 

impedance, Eq* 11 becomes 
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where 


I- * 
mn 


f / 1 


’® 2 <%> ( r > i (¥-) 
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dr 


A related but different expression is obtained for a radially- directed slot. 
The damping coefficient is obtained by solving Eq. (11) for <p t the complex 
frequencies (eigenvalues) which satisfy the equation ; 
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( 14 ) 


where w is the angular frequency of the mode and a Is the damping coefficient 
for tha t mo de. — Conventional root-finding techniques (Newton's method) are 
employed to solve Eq.(l2)but a cavity impedance expression is required. 


The wave motions in the cavity and in the main chamber are coupled through 
the impedance expression, which is an analytical representation of the wave 
motion in the cavity. A constant (uniform) cavity impedance is convenient 
to use but, nevertheless, it is an approximation. Alternately, rather than 
introducing the impedance, continuity of the oscillatory pressure and velocity 
(normal component) can be required as done in the analysis of baffled chambers . 
(Ref. 6). However, this further complicates the analysis and, therefore, the 
cavity impedance has been used. 
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Acoustic Impedance of the Cavity 

Thus, through the use of appropriate impedance expressions, the damping to be 
expected from any form of acoustic resonator may be calculated. However, 
the resonator functions in an environment enormously different from most 
acoustic applications; to be appropriate the impedance expression must per- 
tain to-thls environment. Appropriate impedance expressions are available 
for Helmholtz resonators because extensive work has been done with these 
resonators for use in acoustic liners, 

Must less relevant work has been done with quart erwave resonators and even 
less has been done with generalized variations of these. Hie earliest re- 
ported impedance measurements relating to quart erwave resonators were ob- 
tained during the IH ascent engi ne a coustic cavity study (Ref. 1). Apparently, 
no other data have been published. Recently, some experimental activity has 
been reported as in progress at Princeton (Ref. 7), but no data are yet 
available. In addition, a theoretical analysis was briefly described in 
Ref. 7 to obtain a "time-averaged" acoustic admittance. Because this time- 
averaged admittance was not in a form amenable to ready comparison with the 
Rocketdyne data, no attempt has yet be en made to use it. 



* 

« 
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TWo important experimental observations affect the development of impedance 
expressions for acoustic cavities: (l) the impedance is nonlinear (varies 

with the local oscillatory amplitudes) but, (2) the nonlinearity is princi- 
pally confined to the resistive portion of the impedance. The second obser- 
vation is based on the extensive results available for Helmholtz resonators; 
these results indicate the reactance portion of the impedance is only weakly 
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dependent on amplitude* Furthermore, the nonlinear resistance has been found 
to be s Imply proportional to the amplitude of the oscillatory Mach number in 
the aperture* These two observations were used during the current program 
and during the 111 ascent engine cavity study to formulate expressions for 
quarterwave resonators and later for more general configurations. The fol- 
lowing paragraphs expand upon these remarks. 


For a Helmholtz resonator, the commonly used specific impedance expression 
may be written as (Ref. 8): 



(15) 


where u is the angular resonant frequency of the absorber. £_ is the 
effective length of the aperture and r ^/c. This im p e da nce exhibits 

an amplitude dependence (nonlinearity) at amplitudes of interest for combustion 
stabilization. Fortunately, that nonlinearity is confined, principally, to 
the specific resistance, e^. The linear reactance expression (the imaginary 
term) pertains at high amplitudes with some adjustment in the effective 
aperture length, ? e (which also affects n ). 


In the amplitude range of interest, the linear resistance (effect of viscous 
and thermal-losses) is negligible relative to the nonlinear resistance. Fur- 
ther, in the absence of steady flow through or past the aperture of a Helmholtz 
resonator, the nonlinear resistance is proportional to the amplitude of the 
oscillatory Mach number (peak Mach number) in the aperture. This amplitude 
dependence has been confirmed both experimentally and analytically (Ref. 9). 
This nonlinear resistance, associated with an energy loss, appears to arise 
from oscillatory jet losses at each end of the aperture . That is, an oscilia- 
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tory Jet is formed in the aperture which tends to degrade oscillatory energy 
into turbulence or some other acceptable form. The positive resistance 
implies that osc illat ory energy is removed from the mode of interest; but. 
it is unnecessary for current purposes to know the form of the converted 
energy as long as it is not that of a new Instability mode. 


Thus, the nonlinear impedance of a Helmholtz resonator is composed of a 
nonlinear resistance, due to Jet losses, and a linear reactance, at least 
in form. Because a similar oscillatory Jet should be formed with an acoustic 
cavity (quarterwave resonator) , it is certainly reasonable to expect the 
impedance of the acoustic cavity to exhibit similar impedance character. 
Therefore, the reactance .expansion can be obtained by solving the wave equa- 
tion for the cavity. By adding a nonlinear resistance term, a specific- 
impedance expression is obtained. For a partitioned, straight slot it is 


C - 



(Pc), 

pc 


cot 




( 16 ) 


The specific resistance, e is written as 



(17) 


where Mg is the peak Mach number at the open end of the slot and T is re- 
garded as an empirical coefficient. For Helmholtz resonators r » 1,0, 

roughly; because a Jet can form at both ends of the aperture with a Helmholtz, 
resonator, but only one end with a slot, a value of r~ 0.5 was anticipated 
for the slot. A value of r °* 0.6 was obtained fraa subscale (T-bumer) 
measurements (Ref. 1). Furthermore, the T-bumer data supported the- functional- 
relationship shown in Eq. 16. 
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The resistance expression shown in Eq. 17, with the indicated values of r, 
is appropriate when no steady flow exists either through or past the slot. 
Based on empirical case results from Helmholtz resonators (Ref. 10) the effect 
of steady through-flow is expected to be, approximately 

(Pc)_ A 

9 1 * — r < fl s ♦ V (18) 

where M s is the steady through-flow Mach number. According to currently 
used correlations for acoustic liners, cross-flow is ejected to increase 
the value of r, the increase being related to steady cross-flow Mach number 
(Ref. ll). 


Equation 16 is restricted to partitioned cavities because the wave motion 
in the cavity was assumed to be one dimensional. Related equations could be 
developed for the unpartitioned case, however, this case is of little prac- 
tical interest because such a cavity would probably suffer erosion problems 
due to hot_gas circulation. The partitions in the cavity should be spaced 
at a distance which. is relatively small compared with a wave-length, based 
on the main chamber sound velocity. 


At times, it is of interest to use radially directed acoustic cavities with 
uniform-width partitions, in which case of course, the cavities do not have 
a uniform cross-sectioned area. For this case the reactance expression 
becomes 
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where r-^ and r 2 are the inner and outer radii of the cavity, respectively. 
Equation 19 is difficult to use for practical calculations- because the argu- 
ments of the Bessel functions are complex; however, if the Bessel functions 
are approximated by asymptotic expansions for large arguments, Eq. 19 becomes 


(Pc) 

Xl * - “pc cot k s (r 2 ’ r i> 


( 20 ) 


Equation 20 is equivalent to the reactance expression used in Eq. 16. There- 
fore, Eq. 16 is regarded as suitable for radially directed as well as axially 
directed cavities. 


Eq. 10 was developed with the assu mption tha t the density and sound velocity 
(temperature) were uniform throughout the slot, but different from the main 
chamber values. The same assumption is used in acoustic liner work. However, 
thermocouple data from the full-scale firings show that a temperature gradient 
often, but not always, occurs from the open-to-closed ends of the slot. Less 
severe temperature variations from cavity- to- cavity were also observed. 

These kinds of temperature variations alter the cavity impedance and compli- 
cate Idle analysis. If the cavity temperature is assumed to vary from cavity- 
to-cavity, but to be uniform in each individual cavity, a characteristic 
equation may be obtained from Eq. 11, as before, but the resultant equation . 
contains a double summation rather than the single sum obtained previously 
(Ref. 1). In principle, this case is not different from_the uniform- impedance 
equation; however, it is far more time consuming to solve. 
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If the temperature is assumed to vary from the open to closed ends of the 
cavity, but not from cavity- to- cavity, the u niform -impedance characteristic 
equation. Eg. 12, may be used with a modified Impedanoe expression. However, 
an appropriate impedance expression is more difficult to obtain. The general 
temperature distribution cannot be handled in a simple manner because, for 
steep temperature gradients, the conservation equations cannot be linearised 
in the usual manner to obtain the wave equation. Nonetheless, two limiting 
cases can be handled in a tractable manner: a weak temperature gradient and 

two regions of uniform temperature separated by a discontinuity. In each of 
these cases the wave equation can be used. 

For the discontinuous case, the wave equation may be approximately solved 
employing the long wavelength approximation to give a new reactance expres- 



where allowance has been made for the possibility of a change in cross- ’ 

sectional area as well as temperature at a distance L from the closed end 
of the cavity. The overall slot depth is L + t . . 

4 m 

j \ 

Perhaps the simplest week-gradient- case is. a linear variation in sound 
velocity. The wave equation may be solved to give still another reactance 
expression. 
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where 
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and the sound velocity varies linearly from at the closed end to c fl 
at the open end of the slot. 


Although neither of these temperature distributions accurately represent the 
actual temperat ure var iation , the impedance expressions obtained from them 
should be sufficient to show the qualitative, at least, effect of the temper- 
ature variation. 

Hat mint ed Damning 

Because the specific resistance is amplitude dependent, the linear analysis 
of the combustion chamber described previously is not strictly valid. Never- 
theless, the nonlinearity appears sufficiently weak that a quasilinear inter- 
pretation is satisfactory. Thus, the damping coefficient is calculated with 
the linear equations, but the result depends on the local amplitudes at the 
cavity entrance. Moreover, the amplitudes should be restricted to levels 
less than 20 percent (peak to peak) of chamber pressure because the wave 
equation is not adequate for higher amplitudes (Ref. 3)* These approxi- 
mations are not regarded as highly restrictive but should not be ignored. 

Damping coefficients have been calculated for a variety of cavity, configura* 
tions from the model comprised of foregoing chamber and cavity equations. 

typical results are shown in Fig. 2 for an art ally di rected slot located_at 

the periphery of the injector. These curves were calculated during the 
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earlier investigation, the LM ascent engine acoustic cavity study (Ref. 1). 
Ihe equations apply to either standing or spinning modes of instability; the 
same damping coefficient is obtained in either case. 


For these and subsequent calculations, the specific resistance of the cavity 
has been evaluated -in terms of the local pressure amplitude. This has been 
done in the following manner; 



6 + JX 


( 23 ) 


A ® || 

but & - rM g and let X = X + j x . Ihe amplitudes of the complex 
quantities are obtained from their absolute values; thus 



By rearrangement 


(TM S ) 4 - 2 X"(tA )3 + (x" 2 + x ,2 )(TM ) 2 - (I ±) 2 = o 

a yp 0 


( 25 ) 
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The last equation is solved by a numerical root finding method to determine 
the peak Mach number in the aperture and, therefrom, the specific resistance. 
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FACTORS INVESTIGATED 


At the completion of the previous investigation the model for cavity damping 
was , regarded as largely satisfactory. However, it did not adequately pre- 
dict the optimum slot width, i.e., the cavity width corresponding to maximum 
damping. Full-scale stability rating tests had been made, with cavity con- 
figurations as shown in Fig. 1, from which it was found that 0.3- to 0.4-inch 
wide slots were more effective than 0.2-inch slots. As shown in Fig. 2, the 
model predicted maximum damping with ~ 0.2- inch wide cavities. Consequently, 
much of the analytical effort during the current program has been directed 
toward testing analytical approximations used in the model development and 
removing limitations from the model in the hope of improving the agreement 
between predicted and observed optimum cavity widths and generally strengthen- 
ing the analytical approach. In addition, the influence of a number of 
parameters on stability was investigated analytically. 

Convergence Considerations 

The characteristic equation for an axially directed slot was initially solved 
with an approximate evaluation of the Bessel function integrals. This equation 
was (repeated for convenience): 



(26) 
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where 



Because the slot was generally narrow relative to the variation of the 
Bessel functions for low values of n, approximation of these integrals by 
their integrands multiplied by the slot width appeared reasonable* Thus, 
the following approximation was tried (for all values of n): 



( 28 ) 


A cursory (numerical) examination of the convergence of the corresponding 
approximate characteristic equation indicated that the equation did converge 
and that the approximation was satisfactory* Die results reported in Ref* 1 

were based on that approximation. However, a more thorough investigation of 

convergence was made during the current program; the approximate equation 
was found to exhibit asymptotic character and, indeed, failed to converge 
even with very large numbers of terms, although the term by term variation 
was small* Die approximation was subsequently replaced by an exact integral* 
Comparison of the exact and approximate results indicated a modified approxi- 
mation could be made, i.e. 
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where r g is the mean slot radius and w is the slot width. Effectively, 
this is the approximation used in previous calculations with rl ~ 20. 
Significant errors occur only with large values of n, ~ 100. Nonetheless, 
the exact form of the Qessel function integrals has been used for all sub- 
sequent calculations. 

These considerations led to the conclusion that adequate convergence is 
obtained when d » r w /2w terms are retained in the summation. Note that 
the number of terms required becomes infinite as the slot width, w, goes to 
zero, as should be expected. However., this does not represent a restriction 
because only ~ 20 terms are required at. the optimum width. However, this 
improvement did not significantly change previous results and consequently 
did not remove the disparity in optimum widths. 

End Corrections 

Although unrelated to the optimum width considerations, another factor of 
seme importance is the so-called length end correction that is frequently 
added to the physical length of an aperture to obtain the effective length 
needed for use in the impedance equations. This end correction is usually, 
attributed to a small quantity of gas external to an orifice which is carried 
along with the oscillating gas in the orifice (Ref. 11). With this explana- 
tion of the need for an_end correction, then such a correction should also 

be used in the impedance equations for a quarterwave resonator as well.. 

However, closer examination showed this to be incorrect. 


According to Morse and Ingard (Ref. 4, Pg. 760)1 the correction accounts 
for the flow convergence (and divergence) at the open end of the slot. The 
correction arises because one-dimensional flow equations are used whereas 
locally the flow exhibits multidimensional character. Moreover, the end 
correction is generally calculated from an integral fomu lation similar that 
in use for analysis of the main chamber. Thus, the conclusion was reached 
that the end correction was implicit in the analysis being used and that no 
further correction was necessary for the_damping calculations being per- 
formed. It may be appropriate to include a small adjustment to account for 
nonlinearities in the cavity reactance; however, no such nonlinear correc- 
tion h as been made because no means are available for estimating its magnitude. 


Iterative Approximation 

The accuracy of the variational appro ximat ion, used in the initial form of — 
the model, was verified by comparison with more accurate results obtained 
with a combined variational-iterational method. The iterative method was 
developed from the integral expression for the pressure, i.e., 

p(r) = - Jky |G(r|r s )p(r s )dS Q (30) 

a_ 

where the substitution N* vf> = -Jkyf> has been made. A standard method 

for attempting solution of an integral equation is to insert an initial 

approximation within the integral end perform the integration to obtain an 
improved approximation (Ref. 5, pg.1006). Thus, if 4^jis.the ith approxima- 
tion for the pressure distribution at the slot entrance, the i + 1 approxima- 


t±on_l8 given by 
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Employing this notation the characteristic equation may be simply written as 


/ s (l) (s U) - « (l+1) } 


dS 


( 32 ) 


These equations, for an axially directed slot, were programmed for numerical 
solution. Hie previously used approximate pressure distribution (i.e., the- 
no-slot expression) was used as a_starting point, a zeroth iteration. The 
explicit form of the equations will be shown later under the discussion of 
gain/loss boundary conditions. Examination of the results for several cases 
shows that the original variational approximation was excellent. The error 
with the zeroth iteration in the damping coefficient was less than 0.2 per- 
cent. Furthermore, convergence of the iteration scheme described above was 
found to be excellent. Results from one set of calculations are shown in 
Table I. Hie iterative method implies that the integral equation is solved 
exactly whan ^ i+1 ^ * ^ or £^ i+1 V = 1.00. 

TABLE I . CONVERGENCE OF ITERATION SCHEME* 



Absolute Value 



r = 3.695, 

r = 3.795, 

r = 3.895, 


i 

inches 

inches 

inches 

a, sec" 1 

0 

1.04573 

0.99573 

0.97615 

212.343 

1 

1.00075 

0.99997 

0.99935 

212.684 

2 

0.99991 

1.00003 

1.00006 

212.676 

3 

1.00002 

1.00002 

1.00002 

212.657 

4 

0.99999 

0.99999 

0.99999 

212.672 


*0.2 x 1.75 inch slot, LM ascent engine, p/vp rt » 0 . 06 

F » 0.6, T/T - 4.00. 0 

c s 
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This table Indicates the worst error in the pressure distribution without 
Iteration is less than 5 percent. Further, the error is so small, even with 
one iteration, that it is affected by the accuracy of the root-finding pro- 
cedure. Similar results were obtained for other slot lengths. Clearly, the 
iteration is. not generally needed for damping coefficient calculations. 

A second root of the. characteristic equation was found with the iteration 
method. But this root was found to be trivial because it corresponds to 
4^ = = 0. The procedure converged to essentially this value in 

3 or 4 iterations. 

Hie improved accuracy obtained with the iterational procedure did not 
significantly alter the predicted optimum slot widths. This is shown by the 
calculated damping curve shown in Fig. 3 which were obtained with one itera- 
tion. 

Impedance Averaging 

Another aspect of the calculation which may affect the predicted optimum 
slot width is the manner in which the cavity impedance is averaged. Impe- 
dance averaging is necessary because a uniform impedance approximation has 
been used in the damping analysis. However, because the impedance varies 
with the local pressure amplitude, it cannot be uniform for the tangential 
modes of interest. Previous calculations were made by using a simple root- 
mean-square average pressure amplitude in the impedance calculation. Pre- 
liminary calculations (Ref. 1) indicated this approach was satisfactory. 

The effect of this approximation on the predicted optimum width was-checked 
by comparison with results from a more nearly appropriate average impedance. 
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According to the work of Cantrell (Ref. 3) the impedance should be weighted 
by the square of the local pressure amplitude in the damping calculation. 
Accordingly, a new average slot impedance was calculated from 
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where the circumferential variation of pressure has been taken as 
and the local admittance, y( 0 ) , depends on angular position through this 
pressure variation. The admittance-containing integral was evaluated 
numerically. The resultant average was employed in the characteristic equa- 
tion to calculate the corresponding damping coefficient. Results are shown - 
in Fig. 4 and 5 for two pressure amplitudes. The qualitative character- 
istics of the predicted damping curves are the same as those from the earlier 
simplified calculation. However, .the results indicate an average amplitude 
of ~ 81 percent of the peak value should be used rather than the 70.7 per- 
cent corresponding to an rms average. The optimum slot width is unchangod. 


Reference Amplitude 

Still another modification to the analytical approach was tried* deviously, 
calculations were made with a constant reference pressure amplitude .chosen - 
at the slot, this amplitude being treated_aa a parameter in the-calculat ions . 
A reference is necessary because the cavity impedanee changes with amplitude. 
Because the slot was assumed to be located at the periphery of the injector, 
this amplitude assumption implied a roughly constant pressure (and transverse 
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Figure 5. Predicted Cavity Damping With Integrated Average 
Impedance 
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velocity) at the injector_face and, in turn, relatively constant combustion 
driving# An alternative point of view is to fix the pressure amplitude at 
the nozzle entrance. However, when this alternative was tried, roughly the 
same optimum slot width was obtained but the damping coefficient curves 
were altered somewhat. 


The amplitude of the slot entrance may be related to the amplitude at the 
nozzle entrance through the integral expression for pressure, i.e. 


mn ml r ' 

p.,cos me -ii 
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y' mn m v mir 
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mn . . 


where the coefficients ^ ^ are determined from the iterative calculation 

mn 

and are defined in the next section for the more general gain/loss case. 

This pressure expression was combined with_previously used expressions for 
the cavity impedance and the characteristic equation to calculate the corre- 
sponding damping* ttie results are shown in Fig. 6 . The predicted H«n p^ e 
is altered significantly; but, nevertheless, the optimum slot width still 
appears near 0.2 inches. 
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Gain/Loss End VfrQ J Boundary Conditions 


The next approach used was to include gain/loss admittance- type boundary 
conditions, for both the nozzle and injector ends of the chamber, in formu- 
lating the analytical model. The characteristic equation thus obtained for 
an axially directed slot, is: 
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where 
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In these calculations f the nozzle admittance (y^) and the combustion— end 
admittance (y^) were selected as constant values. Damping calculations were 
also made with this equation and) further) with the reference pressure 
amplitude_chos en at the nozzle entrance. The amplitude at the slot was 
calculated from 
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Premising results were obtained from these calculations) some of which are 
shown in Fig. 7 and 8 . These calculations are in better agreement with 
the observed dependence of stability on slot width than any previous approach. 


The admittance values for the two ends of the chamber have been chosen some- 
what arbitrarily) although the magnitude was chosen based on the so-called 
virtual admittance of a short nozzle. The virtual admittance (Ref • 3) 
includes a convection term, associated with the steady flow, which tends to 
dominate its value; the virtual admittance of the nozzle changes by only 
10Lto-15 percent if the actual admittance is set equal to zero. Further) 
the signs have been chosen so that the injector boundary has a driving effect. 
An alternative would be to replace the injector with a sensitive, time lag 
expression for concentrated combustion) which would be a simple change. The 
injector admittance (virtual) would be simply replaced by 

- joJT 

y i = * M 0 n(1 " e ) (38) 


where n and t are the interaction index and the sensitive time lag) respectively. 
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Parametric Study 


Hie influence of several cavity-related parameters has been indicated in the fore- 
going discussion, i.e., width, depth and, to some extent, amplitude. In addi- 
tion to these, the influence of variations in the resistance coefficient r and 
in the average cavity gas temperature were investigated for a xially directed 
slots. Parametric curves are shown in Fig. 9 and 10 to illustrate the effect 
of varying F . The influence of this parameter was also reported in Ref. 1. 
Gene rall y, the predicted damping increases with r, although for sufficiently 
large values a reversal in this trend is expected. 

The principal effect of variations in the average cavity temperature is one 
of tuning, i.e. the slot depth corresponding to maximum damping shifts roughly 
as-the quarter wavelength varies. However, this effect cannot be eliminated 
by using a nondimens ional slot depth which includes the local sound velocity. 

These e£fects_are illustrated in Fig. 11 through 14. The first of these has 
been reported previously in Ref. 1 but is repeated for clarity. 
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Predicted Cavity Damping: First. 
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Figure 10. Effect of F on Predicted Cavity Damping 
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Radially Directed Cavities 


This gain/loss modification of the model appears to account for the observed 
width effect* However, these calculations of cavity damping have all been 
made with equations pertaining to an axial slot* It was assumed that the 
differences between an axial slot and a radial slot, if each was located 
along the periphery of the injector, wou ld, be small. This assumption was 
based on inspection of the Cantrell equation (Ref. 3 ) and of the integral 
form of the characteristic equation, 

/£y£dS + jkj^y J G(r|r o )y|dS o dS = 0 (39) 

where | is the approximate oscillatory pressure which approaches the actual 
pressure with iterations* This pressure distribution was expected to be 
substantially the same for each slot direction (if the slot is along the 
periphery of the injector) • Therefore, the slot damping was expected to be 
substantially the same. Consequently, the axial formulation was used because 
it is somewhat more tractable from a numerical standpoint. This assumption 
has been confirmed* contrary to initial reports. 

When the radial slot formulation was programmed so that the influence of the 
axial position of a radial slot could be investigated, significant (but 
erroneous) differences were initially found and reported in the monthly 
reports and in Ref. 12. Subsequently, an error was found in the computer 
program; the cavity admittance was incorrectly multiplied by the slot width. 
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Corrected results show the calculated damping to be nearly the same for both 
axially and radially directed slots located near the periphery of the 
injector. 


For a radially directed slot and with no. end-wall gains or losses, the itera- 
tive fom of the characteristic-equation is 


1 4ru- ¥ ¥ t ■ 0 < ko > 

where 



and 



1.0 q = v 

0 q i v 


The Bessel function ratio is being evaluated from the following expansion 


Jg( y y) = r 5* 2Q ^ 

w ~ ’- 2 




( 41 ) 


This expansion w rs obtained in Rsi. 8. 
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The foregoing equations have been programmed, for numerical solution as done 
previously for an axial slot, with the reference amplitude taken at the slot. 
Calculated damping £or_twQ different reference amplitudes are shown in 
Fig. 15 and 16. For these calculations, the upstream edge of the slot was 
assumed to be adjacent to the injector. Comparison of these results with 
the corresponding damping curves lor axially directed slots (adjacent to 
the chamber wall) show only minimal differences. The calculated effect of 
moving a radially directed slot downstream is shown in Fig. 17- The param- 
eter in these calculations was the axial location of midwidth of the slot, 
which is indicated on the figure. The results indicate that the maximum pre- 
dicted damping is only weakly dependent on axial position but this maximum 
occurs at different slot lengths (tuning). The latter effect is somewhat 
surprising; it must occur because of the oscillatory interaction between the 
cavity and main chamber. 

The results indicate an optimum slot width near 0.2 inches irrespective 
of whether the slot is radially or axially directed. Thus, the predicted 
optimum slot width, without the inclusion of gain/loss boundary conditions, 
does not agree with the full scale stability results. Some indication of 
the reason for this disparity may be obtained from the erroneous results 
mentioned earlier, which appeared to agree quite well with observed engine 

stability. These calculations were made with the slot admittance 

weighted (multiplied by) the slot width, i.e. y - wy/w 0 where w Q la. of 
unit width. The corresponding results are shown in Fig. l8 and they agree _ 
quite well with observed engine stability. 
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Figure 18. Predicted Slot Damping With Area-Weighted Admittance 
for Radially Directed Slot Adjacent to Injector 
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The observation that the predicted damping is more closely in agreement with 
measured stability when the admittance is weighted by the slot width, suggests 
that the optimum slot width is underpredicted because the impedance being used 
is too low. The absolute value of the impedance may be increased by in- 
creasing either the amplitude or the coefficient r, because the specific 
resistance increases with either of these. Therefore, several calculations 
were made with substantially larger values of r to investigate its effect 
on the predicted optimum slot width. Results from these calculations for 
values of 1.1, 1.5 and 3.0 are shown in Fig. 19 through 21 . Clearly, the 
optimum .slot, .width shifts in the proper direction wit h in creasing r • In- 
deed, the predicted curves for r = 1.5 and 3.0 are in good qualitative 
agreement with the stability data. Although these values of r are in 
excess of those found in the T-bumer tests, the results can be explained as 
follows. 

First, it may be observed that the effects of increasing r or of pressure 
amplitude at the slot are interchangeable. These parameters appear in the 
analysis in a combined form, i.e. p^/yp o . Thus, the agreement indicated 
above between predicted and observed slot widths can also be obtained by 
increasing the value of the amplitude at the slot used in the calculations. In 
addition, as indicated in Appendix A, the effect of steady cross flow in the 
chamber (for radial data) is also to increase r. Presumably, secondary 
flows (recirculation) in the engine could have a similar effect. However, 
based on the most recently reported work by Garrison on the effect of cross 
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flows with Helmholtz resonators (Ref . 11), the increase in r probably would 
not exceed 50 percent, i.e., r =0.9. Thus, a better explanation appears 
to be the average amplitude used in the calculations. 

Previous calculations have been restricted generally to chamber amplitudes 

less than 20 percent of chamber pressure, peak-to-peak, but only becau se the 
wave equation loses validity at higher amplitudes. The fact that the best 
agreement is obtained with larger assumed amplitudes indicates that these _ 
should be used in design calculations > but it suggests even better predictions 
of cavity damping could be obtained from a fully nonlinear model. 

TEMPERATURE DISTRIBUTION WITHIN THE SLOT 

The acoustic properties of the gases within an acoustic cavity are of great 
importance to the cavity damping processes. To facilitate the accurate pre- 
diction of the cavity temperature distribution, a semi. empirical analytical 
model was developed from the slot temperature data obtained during this pro- 
gram. The analysis was undertaken, in part, to attempt adaptation of some 
existing model, for heat transfer with flow over a slot, in the hope of predict- 
ing the cavity temperature distribution for any propellant combination. 

Unfortunately, no suitable model was found in the literature. The model 

obtained does not allow accurate prediction of the temperature distribution 
for propellant and/or injector configurations different from. the_one used 
during the hot firing program. Because the analysis is of limited utility 
and, further, the discussion is relatively lengthy, it is included as an 
appendix rather than here (Appendix B). 
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UNCONVENTIONAL CAVITIES 


Considerable-attention has been given during the current program to what have 
been loosely termed as unconventional acoustic cavities. These are acoustic 
resonators of more general shape than either simple Helmholtz or quarterwave 
resonators. These more general resonator shapes are of Interest because 
frequently Pimple, well-defined acoustic resonators of either the Helmholtz 
or quarterwave type are not practical In an engine due to hardware or spatial 
limitations. The h ar d ware chara cteristics often suggest resonator configura- 
tions which do not clearly fit into either category. Moreover, frequently 
neither of these resonator types represents the optimum configuration in 
terms of rnqyimum damping and/or bandwidth attainable within a specified volume. 
Consequently, work has been done to extend the range of resonator configura- 
tions which can be described analytically. Further, some experimental-work 
has been done to characterize these unconventional resonator configurations. 
Some of the configurations under consideration are shown in Fig. 22. Hie 
intermediate resonator is so termed because it is intermediate between a 
Helmholtz resonator and a quarterwave resonator. The intermediate, "L“- 
shaped, Helmholtz and quarter-wave resonators can all be considered special 
cases of the generalized resonator. 

Acoustic Impedance of Unconventional Cavities 

Acoustic impedance expressions have been developed for each of the configura^- 
tions Shown in Fig. 22. Ihe same analytical approach has been used as was 
used for the quarterwave slot, i.e., combine a nonlinear empirically-defined 
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resistance with a linear reactance expression. Furthermore, the same form 
of the specific resistance has been used, 

o « rA s (42) 

where the coefficient P must be measured for these cavity configurations. 

An ap proximate reactance expression can be obtained for the intermediate 
case through use of the long wavelength approximation wherewith the pressure 
and volumetric flow at-the. ...area„change in the resonator are equated. The 
resultant expression is 


„ _ uiL w 
<*>. C0S c s 


(pc) S 

S sin — sin 
' pc) c S s c c c s 


PC 


“* tin 5i „ cot 


( 43 ) 


where allowance has been made for the possibility of a change in temperature 
as well as cross-sectional area at a distance L from the closed end of the 
cavity. The overall cavity depth is L + fi . This expression limits to 
the quarterwave expression when (pc) c S c — (pc) fi s 3 ... It limits to the Helmholtz 
resonator expression when the cavity dimensions are small compared to a wave- 
length} coC/c << 1, coL/c « 1 so that small-argument approximations may be 
used for the trigonometric functions. Moreover, this expression indicates 
that a quarterwave resonator occupies a minimum volume relative to inter- 
mediate or Helmholtz resonators with the same resonant frequency* This-result 
is shown in Fig. 23 as calculated from this equation, Eq. 43. 
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A more general reactance expression wa s devel oped for the generalized 
resonator which includes the Above cases as specif?! cases. A resonator 
shape of the following form was analyzed: 



Through application of an analysis (R ef. 6) de veloped to describe the. wave 
motion in baffled chambers, the following reactance expression was obtained, 
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This equation corresponds to a zeroth iteration. A more nearly exact re- 
actance value can be obtained by solving the following implicit equation 

2 ^ [V“ ( V * V] a q 1+l) ] * 0 (**5) 

where tan = Jk/k q X 
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0 q / 0 


The additional complexity of the latter method prevents its u se -i n .-mos t 

cases. 

T-burner tests which will be described subsequently, have been made with 
both “L"-shaped and intermediate-type resonators to measure the nonlinear 
resistance of these resonators. The results were correlated with the same 
form of relationship for resistance as that used for straight slots. The 
experimental results are scattered but indicate larger values of the empiri- 
cal coefficient, r , than that obtained for quarterwave or Helmholtz 
resonators. The experimentally determined values for r were 1.8 and 4 for 
intemediate and "L"-shaped resonators, respectively. 

Resonant Frequency Calculations 

The resonant frequencies of these resonators indicate roughly the frequency 
ranges over which they will provide damping. The resonant frequency, (actually 
the undamped resonant frequency) is the frequency at which-the cavity, re- 
actance is zero* Therefore, it may be readily calculated from the foregoing 
expressions. 
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Resonant frequencies for two unconventional configurations are shown in 
Table 2 , The lowest two resonant frequencies are shown for the "L" shaped 
resonator because these roughly correspond to quarterwave and Helmholtz 
resonances. This result suggests the possibility of obtaining broad band- 
width by tailoring the proximity of these two resonant frequencies or 
damping two modes of instability with a single resonator by tuning the two 
resonant frequencies to damp these modes. 


Damping Calculations 

The cavity damping to be expected from unconventional cavities may be cal- 
culated as before by substituting these impedance expressions.- Damping cal- 
culations have been done for resonators of the "L"-shaped and generalized 
types. The resultant damping coefficient curves are largely similar to those 

obtained for straight slots. However, greater bandwidth is indicated in - 

some cases. Results from some of these calculations will be described sub- 
sequently in connection with the cavity design for the full scale tests 
with unconventional cavities. 

OTHER FACTORS 

The contractual work statement for this program specifies a number of- cavity 
and stability related factors that were to be investigated^ Hot nil .of 
these could be investigated in detail within the scope of the program. Other 
portions of this report describe the detailed investigations that were per- 
formed. In this section, those factors whieh could not be Investigated in 
detail will be discussed. 
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Steady Flow Effects 

Generally the effects of steady flow have been ignored in the foregoing 
analyses, the assumption being Implied that these effects may be included in 
the composite (negative) damping coefficient provided by the combustion 
system. To investigate the importance of this omission, a variation of the 
model was partially developed which explicitly included a uniform steady 
flow in the chamber. Unfortunately, this modified analysis could not be 
completed within the scope of this program because the steady flow considerably 
complicates the analysis. Some of the equations were developed but no numerical 
results were obtained. 

For the uniform steady-flow case, an inhomogeneous wave equation must be 
solved, it being 

V 2 p + k 2 p - 2jk M o JLfi + MfJ 2 -^§ (46) 

dZ 

The boundary conditions are 

y p H* 

M • VP = - JkyjS M o N (1.7) 

Where M q is the steady flow Mach number and y is the local specific acoustic 
admittance. The inhomogeneous wave equation can be readily converted to an 
integral form but solution of that integral equation is complicated by the 
boundary equation. A major Simplification results i f a ve locity potential 
is introduced because the boundary condition snters-in a more tractable form. 
The velocity potential satisfies the same differential equation as that shown 
above for the pressure) however, the bounda ry condition may be directly 
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written in terms of the admittance and pressure. Converse ly f in the boundary 
condition shown above for the oscillatory pressure, the oscillatory velocity 
at the boundary is not eliminated by introducing the admittance. 

Because of these kinds of difficulties, the steady flow model was not completed. 

Interaction Between Baffles and Acoustic Cavities 

No analytical model is currently available for the analysis of cavity effects 
in baffled chambers. However, such a model can be directly developed by 
combining the cavity damping model described previously with the Rocketdyne 
model for the analysis of baffled chambers (Ref. 6). Ihe latter model has 
only recently been extended to include cylindrical chambers. Although these 
models can be conceptually combined in a straightforward manner, the mechanics 
of doing so is complicated. Moreover, the resultant computer program would 
probably require long computing times. Therefore, the combined model has not 
been developed. 

Nonetheless, calculated results from the individual models indicates no reason 
to expect an incompatibility between baffles and cavities. Certainly none 
exists when the cavities are designed to damp instability modes which are 
compatible with the baffle, e.g., the third tangential or first radial modea 
occurring with three equally spaced radial baffles* The flight configuration 
LM ascent engine employs this combination with complete success. In this case, 
because the mode is compatible with the baffle, it is unaffected by the baffle 
and the existing cavity model can be used directly* 


The alternate case^in which the mode is affected by the baffle, would require 
a combined model for analysis. However, individual-analysis of the acoustic 
behavior of baffled chambers and measurements with acoustic models Show that 
the normal modes of the chamber are distorted by the baffles such that the 
pressure amplitude is increased at the baffled end relative to the opposite 
end of the chamber. This effect of baffles is likely to augment the damping 
produced by a -cavity located near the injector face. _ 

Thus, although specific analysis has not been done, the available information 
indicates baffles and acoustic cavities can be effectively used together and 
that baffles may complement the cavity. However, some caution should be used 
because the stabilization mechanism of baffles has not been established; the 
cavities may tend to degrade the baffle proc esses if these involve nonlinear 
damping effectsu. 

Effects of Varying Engine Design 

Among the other factors of interest are those which affect the applicability 
ofLacoustic cavities to other engines. Specifically what are the capabilities 
and limitations of cavities when chamber diameter, thrust level, energy release 
rate, etc. are changed? Generally, these factors cannot, be evaluated thoroughly 
without either hot firing the configuration or use of an instability model. 
However, in the context of the currently, developed cavity damping, model these 
factors may be regarded as simply affecting the driving processes {which -are 
uns pecified in tha model), except for the chamber dimensions aad_gas properties 
which appear in the model. Thus , the damping model is insensitive to these 
changes and applies equally well to all of them. 
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Other evidence which suggests that cavities or acoustic-absorbers in general 
may be applied to any configuration is simply the wide range of engine con- 
figurations which have already been stabilized with absorbers (see Appendix A). 
Nonetheless, there are evident limiting factors. First, as the 
chamber diameter, is increased the volume requirements for tuning a cavity 
become large. Moreover, the cavity damping coefficient decreases_with__ 
frequency or as the inverse of chamber diameter while, qualitative, the driving 
processes not to decrease as rapidly. These suggest greater difficulty in 
stabilizing a large engine than a small one, as might be expected. 


Another .specif ic factor of. interest is injector shape, i.e., flat or concave. 

This factor may affect the energy release characteristics, the mode shape to a 

limited extent and the proximity of a cavity to the injector "face" or energy 

release patt ern. None of these are likely to affect cavity performance 

other than to change damping requirements and/or tuning (cavity sound velocity 

and density). Propellant conditioning temperature should have a similar effect. 

Although each of the foregoing factors is of importance and interest, they 
were regarded as of lesser significance and importance than those factors that 
were investigated thoroughly. 

t 

ASSESSMENT OF CURRENT ANALYTICAL STATUS * 

During this program the analytical approach for cavity design has-been 

significantly improved. Many weaknesses have been eliminated but seme still 

exist. When the program began the principal concern relative to the analytical ' 



the slot width or- open area corresponding to max i mu m damping or stability. 
Thus, much of the analytical effort was directed toward this deficiency. 



1 


f 


The accuracy of the damping calculation has been significantly Improved, i.e., . 
the equations which comprise the model are solved more accurately now. The 
convergence characteristics of the series expressions have been examined and 
approximate criteria for convergence have been developed. The uncertainty 
relative to aperture end corrections has been removed. The average cavity 
Impedance has been clarified. A model for radially directed cavities has 
been developed. And, an iterative method has been developed to eliminate 
uncertainty relative to the variational approximation and to allow the 
approximation to be improved to .any desired accuracy. These Improvements 
provide assurance that the equations are accurately solved. 

« 

The analysis has shown that the observed optimum width, which is only weakly 
defined by the full scale data, can be predicted if the cavity impedance 
parameter rp/vp o is made large enough. This constitutes a satisfactory 
explanation of observed behavior. However, if reasonable values for r are 
assumed then the corresponding pressure amplitude is too large to be accurately 
described by linear acoustic equations. Nonetheless, the error introduced by _ 
this nonlinearity is probably not as large as other uncertainties. Moreover, 
it can probably be lumped Into the undefined damping requirements. An alter- 
nate explanation was obtained by introducing gain/loss boundary conditions_at 
both ends of the chamber; however, this explanation is less satisfactory 
because appropriate values for the end boundary conditions are not well defined. 
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In addition, significant progress has been made concerning unconventional 
cavity configurations. A satisfactory model has been developed for the 
acoustic impedance of a rather general cavity configuration. Moreover, it 
was demonstrated that a quarterwave resonator occupies a mi nium volume 
relative to Helmholtz or intermediate-resonators for a fixed resonant 
frequency. 

The analysis of cavity temperature distributions was not entirely satisfactory 
because the resultant model does not allow prediction of cavity temperatures 
in new engine environments. 

Several weaknesses are evident in the analytical approach for predicting 
cavity damping. One is the approximate separation of contributory processes 
through use of the damping coefficient. This approximation loses its validity 
when the contributory processes interact significantly. Moreover, even if 
this approximation is valid, the inability to predict the damping requirement, 
or net driving contribution, is restrictive. In addition, the quasilinear 
approximation has been used to aid tractability, but this must introduce some 
error. The effects of steady through flow have not yet been included; they 
can and should be added. Finally, the acoustic cavity may alter the wave 
motion in the chamber so that the oscillatory driving is reduced j this effect 
cannot be predicted without a model for the combustion. 

Nonetheless, the analytical model for cavity damping has been proved useful. 

In spite of its limit ations, it is believed to be the most .accurate and re- 
liable method available for the design of any kind of acoustic absorber. 


SUBSCALE TESTS 


Two series of subscale tests were made during this program. The first of 
these was made tc evaluate experimentally the potential for increasing the 
frequency range (bandwidth) over which a cavity system is effective through 
the use of multiply-tuned cavity combinations. The purpose of the second test 
series was to measure the specific resistance of unconventional cavity con- 
figurations. Generally, the techniques used are unchanged from those used 
during the LM ascent engine acoustic cavity study (Ref, l). 

For these tests, a combustion-driven acoustic oscillator was used to generate 
the necessary high-amplitude environment. This oscillator was a variation of 
the T-burner which is extensively used in studies of solid-propellant combus- 
tion instability. An acoustic. cavity or a combination of cavities was mounted 
within the T-bumer and the pressure response of the cavity was monitored 
with_high-frequency pressure transducers. The T-bumer provides an oscillatory 
environment for cavity evaluation which more closely simulates the engine 
conditions under which the cavity must operate than a cold-flow impedance 
device. 

The bandwidth tests were made over a range of frequencies with several combina- 
tions of cavities. The damping contributed by each absorber combination was 
inferred from the observed limiting pressure amplitudes. 

For the impedance determination,, pairs of oscillatory pressure measurements 
are made from which the impedance was inferred. The technique m ay be simply 
illustrated for a straight quarterwave resonator, as fellows. 
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The specific acoustic impe dan ce of a cavity is 


f = 



( 48 ) 


For the fundamental frequency component, the oscillatory velocity at the open 
end of the slot, u s , may be related to the pressure at the closed end by 


u c 



These may be combined to give 
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" 0 » = ) s 


£ sin 

c s 


(49) 


(50) 


If the. slot impedance expression, Eq. 16, is introduced, Eq, 50 becomes 

P, 


r ° , ,, wfio 

S j i sin — - + cos ® 


(5D 


Further, the amplitude of the oscillatory Mach number at the slot entrance is 


A 




(52) 


from Eq. 49. % measuring the amplitude ratio and phase angle between the 

oscillatory pressures at the open and closed ends of the slot, the slot 
impedance can be determined from these equations. The same approach was u sed 
for the unconventional cavities. The T-bumer experiment was set up so that 
the necessary pressure measurements could be made. 
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T-BURNER SETUP 


The T-buroer used for these measurements is the same as that used previously] 
it is shown schematically in Fig. 24 . The interior of the motor was cylindri- 
cal, closed on both ends, with a nozzle in the side near the midpoint. An 
annular ring of solid propellan t was burn ed in one end, as shown in Fig. 24,. 
and the acoustic slot or slots to be tested were mounted in the opposite end 
of the motor. With the arrangement shown, the propellant combustion tends to 
spontaneously_excite a first-longitudinal mode of oscillation, the degree of 
excitation being dependent on the propellant used. 

The T-bumer was comprised of interchangeable 1.5-inch ID lengths of mild 
steel tubing, a nozzle assembly, and two end assemblies. T-bumer lengths 
(interior) of 5*0, 9*0 and 13*0 inches were used. These lengths, of course, 
determine the oscillatory frequency experienced during, a test (nominally, 

3400, 1750, and 1100 Hz). One end assembly contained a simple igniter; the 

other contained the acoustic cavity to be tested. The cavities were inters 

changeable so that a wide range of cavity dimensions could be tested* The 
cavities were made of steel to assure dimensional stability and avoid degrada- 
tion with usage. 

The composition of the solid propellant used during the study was: 78 weight- 
percent ammonium perchlorate (AP), 2-percent copper chromite burning-rate 
catalyst, and 20-percent carboxy-berminated polybutadiene binder* A bimodal 
blend of coarse (-200 micron) and fine (- 20 micron) AP in the preportion of 
1:3 was used to promote castability* The propellant was cast, into phenolic 
sleeves and cured. The easting was then machined to length and Internal 
diameter on a lathe. 
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SLOT 







The propellant, was ignited by a s m a ll amount of propellant sawdust, which, 
in turn, was ignited by a short piece of electrically! initiated exploding 
bridgewire (Pyrofuze). Because an annular ring of propellant was used, the 
burning surface area increased as the surface regressed; therefore, chamber 
pressure increased. The nozzle was sized so that the pressure passed through 
120 psia about midway through the test. This ramp in static pressure (from 
~70 to -200 psi) is too slow to affect the high-frequency data. An annular 
ring of prope llant , rather than an end-mounted disk, was used to increase the 
burning area and, thereby, the oscillatory driving. 



At a chamber pressure of 120 psia, typically peak— to-peak amplitudes of 60 psi 
without a slot and 20 psi with a slot were obtained. Oscillatory pressures 
were measured with Kistler (Model 603A) high-frequency transducers. These were 
used because of their stable sensitivity and very high resonant frequency 
(-400,000 Hz) which are required for accurate amplitude and phase measurements. 
The oscillatory-pressure data were recorded directly on an oscillograph and 
on magnetic tape. An oscilloscope was used for "quick look" purposes. 

Considerable care was exercised with the data recording system to assure that 
corrections could be made for disparities in the gain and phase characteristics 
of the two recording channels. Before each test, a reference sig n al, with a_ 
frequency near that expected for the data, was simultaneously introduced into 
both recording channels at the output of the Kistler charge amplifiers. This 
recorded signal was subsequently used to make amplitude and phase corrections — 
to the data; these corrections were minimized by electronic adjustments*— 
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Amplitude measurements ware made dirootly from the oscillograph records. An 
electronic phase meter (Ad Yu, model A 05) was used to make phase measurements 
from the tape-recorded data. Because only the ampHtude and phase of the 
fundamental frequenc y comp onents of the data (rather than the composite signal) 
were needed, the data were electronically filtered for each kind of measure- 
ment. Allowance was made for the phase shifts in the filters. Hie electrical 
output from the phase meter was recorded on an oscilloscope. 


BANDWIDTH TEST SERIES 

Wide bandwidth is a highly desirable feature of any acoustic absorber design. 

If the absorber is effective over a wide frequency range, the possibility of 
eliminating one mode of instability only to encounter a new mode may be pre- 
vented. Further, a wide band absorber is less sensitive to the poorly known 
environmental conditions. Therefore, means for achieving mAviimim bandwidth 
were explored. 

The only evident means of expanding the bandwidth of conventional acoustic 
cavities (straight quart erwave resonators) is to employ parallel resonators 
with different resonant frequencies. However, these absorbers are not inde- 
pendent of each other and the interaction may destroy the overall effectiveness. 
That is, the combination may, under seme conditions, be less effective than a 
single cavity, in terms of maximum damping and actual bandwidth above a 
specific damping level. 


These effects were investigated in a series of subscale tests employing the 
T-bumer. 
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Bandwidth Configurations 


The damping provided by (l) a single cavity, (2) by two parallel cavities and 
(3) by three parallel cavities (the combination cavities having different 
tuning) was measured over a range of frequencies. Each of these cavity sets 
was mounted in one end of the T-burner and tested with three different 
T-burner lengths (frequencies). 

The T-burner hardware is the same as that used previously during the LM 
ascent engine cavity program, except that an improved cavity test section 
was used. This test section is shown in Fig. 25. The principal change is 
in the material used for the_cavity. High temperature plastic was used 
previously, for expediency, whereas a laminated (braced) steel cavity was 
used for these tests, thsb chang e was made to prevent dimensional and flow 
variations with usage and, hopefully, improve the repeatability of the test 
results. In addition, several thermocouples were introduced to improve 
definition of gas conditions (environment). 


An analysis of the damping to be expected from each combination of cavities 
was used to aid selection of the configurations to be tested. A one-dimensional 
analysis was used with the cavity admittance included on an area- weighted basis. 
Two analytical models were formulated by this approach. One mode l co nsidered 
the temperature in the T-burner to be uniform and the other included a step 
temperature discontinuity (T-burner studies indicate such a temperature dis- 
tribution is appropriate). I^Pical analytical results are shown in Figure 26. 
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these calculations show that the cavity damping coefficient is a nearly linear, 
function of pressure amplitude in the T-bumer, for each T-burner length 
(frequency), the results inply that the pressure amplitude In the T-bumer 
will Increase until the cavity damping balances the net driving influence 
provided by the T-bumer. In terms of a negative damping coefficient, this 
driving Influence can be estimated from T-bumer test data with the cavity 
removed. 

Figure 2 6 also shows a comparison between the analytical results an d experi- 
mental results obtained during the previous acoustic cavity program. 3he 
dashed curve for the two-temperature case was calculated by assuming the 
high-temperature section occupied two-thirds of the T-bumer. 

Based on thesa_results, the matrix of cavity configurations shown in Table 3 
was selected, fabricated and tested. 

Bandwidth Test Results 

The assumption used to analyze the bandwidth test data was that, at the 
lim i t ing pressure amplitude, the damping provided by a cavity system was 
equal to the acoustic energy. provided by the combustion processes. The 
limiting pressure amplitudes were measured from the data and the corresponding 
damping coefficients were approximated from the growth-constant data obtained 
previously in blank tests (no cavity) for the three T-bumer lengths. The 
results are shown in Figures 27 through 29 along with analytical results from 
the one dimensional model. 
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TEST MATRIX FOR BANDWIDTH TEST SERIES 












Nondimensional Amplitude, {>/ Y p 0 
Figure 28* T-BuMer Damping Result a, Double_Slot — 






Damping Coefficient 



Nondiraensional Amplitude, y p Q 


Figure 29* T-Bumer Damping Result a, Triple Slot 








The agreement is somewhat less than satisfactory, but the. disagreement is very 
likely associated with the experimental results and not the analytical model. 

The calculations are based on an assumed, and fixed, temperature distribution 
which is uniform throughout the T-bumer. The temperature data, shown in 
Fig, 30, indicates a significant temperature variation} these data were obtained 
from an iron- cons tantan thermocouple mounted near the slot-containing end of 
the T-bumer. 

Generally, the T-bumer slot damping data are believed compatible with the 
calculations, but the variability in slot damping that can occur under irregular 
operating conditions is evident. 

TEST SERIES WITH UNCONVENTIONAL CAVITIES 

Previously, primary attention has been given to conventional straight quarter- 
wave resonators or to Helmholtz resonators. However, often it is difficult to 
find sufficient space in typically designed engines for a straight slot, 
especially if it is-long. Thus, it is desirable to consider some alternatives 
which allow greater design flexibility, A series of tests were made to measure 
the acoustic impedance of severed, practical variations of the acoustic cavity, 

Hie intermediate and "L" shaped unconventional cavities have been described in 
the Analytical Studies section of this report. Six cavity configurations of 
each of these kinds were selected, fabricated and tested in the T-bumer. 

These configurations were selected largely on a basis of their resonant fre- 
quencies, which were chosen to have proximity to the T-bumer frequencies in a 
range that have given the best results in earlier tests with straight slots. 
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Temperature (°F) 



During this test series, 28 T-bumer tests were made to evaluate the impedance 
characteristics of intermediate and of ir L "-shaped cavities. The completed 
test matrix, showing the configurations tested, ia shown on Table 4 . A 
typical configuration is shown in Pig. 31. For each test, the chamber pressure 
and cavity response were measured with pressure transducers and a thermocouple. 
These measurements were then used to determine the acoustic impedance of the 
unconventional cavities. 

Unconventional Cavity Test Results 

The data from the 28 T-bumer tests were analyzed in a manner similar to that 
outlined earlier in this section of the report for quarterwave resonators. 
Generally, this analysis assumes the cavity reactance to be known; the specific 
resistance is then calculated from the data through use of the reactance. ex- 
pression. The accuracy of the reactance expression was qualitatively checked 
by comparison of measured and predicted amplitudes and phase angles. 

For analyzing the data, a generalized impedance expression was formulated to 
relate the cavity-response measurements to the specific acoustic resistance 
of the cavity. The expression is 

d s = - tan 4/ tan 4^ (53) 

where 0 8 is the specific resistance, 4* is the observed phase angle, between— 
the oscillatory pressures at the closed and open ends of the cavity, and 4* 0 
is defined by the following equations 
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The terms in this expression have been described previously in the Analytical 
Studies -section • 


In addition, 

r p cos + = __li (55) 

{ „ 1 e„ k cos m7r A 

sin k£ + cos k£ tan 'F fT'— 221 

0 W k sin k L 
m m 

where k = w/c ft and is the measured ratio of the oscillatory-pressure 
amplitudes observed at the open and closed ends of the cavity. 


Equations 54 and 55 were solved simultaneously for 4* 0 and k corresponding 
to the measured 4* and rp. The specific acoustic resistance was then calcu- 
lated from Eq. 53 . The oscillatory Mach nuuaber (M ft ) at the open end of the 
cavity can be calculated from the following relationship. 
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A linear correlation between the acoustic resistance ( 6 a ) and the peak Mach 
number (M*) was sought, i.e*, 0 a = (a straight line passing through the 
origin). The corresponding calculated results for the "L "-shaped and interme- 
diate cavities are shown in Figures 32 and 33. The scatter in the calculated 
results is not unexpected. Obviously, many other correlating equations could 
be used; however, the form chosen is compatible with available data and models. 
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Figure 32. Empirically Determined Specific Acoustic 
Resistance for Intermediate Cavities 








Therefore, only this equation was used* The l' value for an intermediate 
cavity, based on those data, is approximately 1.8. For a "L"-shaped cavity, 
T is between 3 «0 and 5»0. 


Cavity temperatures were measured with an iron—constantan thermocouple during 
each test* Figure 34 shows the cavity-tenperature results plotted versus the — 
Tr-burner length. 


In addition, qualitative comparisons were made between th e me a sur ed and cal- 
culated values of the amplitude ratios and phase angles,. The calculations 
were made from the equations shown above and with assumed values of the 
parameter: rPj/7P 0 • Results from these comparisons are shown in Fig. 35 
and 3 6 . The results are plotted as amplitude ratios (l/r p ) and phase 

angles (*) versus co/co o with rv ± hl9 0 as the parameter. A single value 
of w 0 was selected by plotting which best fit all the experimental data. 

The agreement between the experimental results and the analytical predictions is 
not very satisfactory but no other model is available for analyzing unconventional 
cavities. Therefore, these expressions were used. 


LIQUID PROPELLANT T-BURNER 

Early in the program the possibility of converting the T-bumer from solid 
propellants to liquid propellants was briefly explored. Several tests were 
made with a single triplet infection element located in one end of the TV burner. 
The propellants were Hrfi^/50% UDMH. Results from these tests were 

somewhat inconclusive because of measurement difficulties, but they Indicate 
the T-burtier was not unstable. No attempt was made to pulse the motor and. 
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T-Burner Length (in.) 

Figure 34. T-Burner Temperature Data From Unconventional 
Cavity Test Series 
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Figure 38. Measured Amplitude Ratio From an "I," -Shaped Cavity 
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thereby, initiate an instability. Because of these discouraging results, no 
further attempts wore made at this conversion. The conversion was tried 
because the use of liquid propellants offered some advantages in cost and ease, 
of operation. 


ASSESSMENT OF T-BURNER RESULTS 

The results from the T-burner tests were somewhat less than satisfactory, 
principally because the data were too badly scattered to draw firm conclusions. 
This unsatisfactory character to the results does not imply that the technique 

is invalid or not useful j rather, it probably implies that greater care must 

be exercised in performing the tests along with seme redesign of the hardware. 

The data scatter is very likely due to nonuniformities and run-to-run varia- 
bility in the burning surface area, which may be caused by inadequate ignition 
or inhibiting along the radial surfaces of the propellant. Such problems are 
not unus ual with T— burners but were expected to be less significant in this 
application. These problems can be eliminated; a more positive and uniform 
ignition method can be used and the propellant can be bonded in place or 
cast into a sample holder. Consequently, significant improvements in the 
quality of the T-burner results appears possible but probably at greater expense. 

In spite of these difficulties, some conclusions may be drawn. _ The bandwidth 
test data are believed generally compatible with the analytical predictions. 

Thus, the model is believed adequate. The results from the unconventional 
cavity tests are more nearly satisfactory. Scatter in the specific resistance 
results is not excessive but the slope of the resultant. correlation is higher 
than appears reasonable, i.e. , the value of f seems too large. The comparison 
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FULL SCALE STABILITY TESTING 



I 


Full scale stability rating tests have been made to evaluate the effects of 
several cavity-related parameters on stability and, conversely, to evaluate 
the influence of other parameters on stabilization with_acoustic cavities. 

Four-different series of tests were performed; these encompassed 125 hot 
firings . 

The first series of tests (Phase I - Full Scale Test Series) was made to 
measure the influence of cavity (l) width, (2) axial position and (3) multi- 
plicity on engine dynamic stability. The second series (Off-Nominal Test 
Series) was made to measure the infliienca _of engine operating conditions, 
i.e., mixture ratio and chamber pressure on cavity stabilization. The third-, 
test series (Film Cooling Test Series) was made to measure the influence of 

l 

film. coolant flow rata on cavity stabilization. Finally, the fourth series ^ 

of tests (Unconventional Cavity Test Series) was made to evaulate the effec- 

1 

tiveness (ability to promote stability) of unconventional acoustic cavities. j 

* ) 

In addition, detailed cavity temperature measurements were made on all tests. \ 

•i. 

1 

; 

The hardware used for these tests was similar to the LM- Ascent engine but with an 
unbaffled injector being used. However, the propellant valves and ducts were not 
similar. The injector is the same as that used previously during the LM Ascent 
Engine Acoustic Cavity Study (Ref. 1); it includes 214 unlike doublet— (pairs ) 
elements and 30 equally spaced, axially directed, film coolant orifices. The in- 
jector was designed for nominal operation with HgQ^/ 50# NgH^ - 50£_UDMH propellants 
at a mixture ratio of 1 . 60 and a chamber pressure of 122 psia . A solid- walled steel 
thrust chamber was used for these tests; it was designed to accommodate replaceable 
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rings which contained the acoustic cavities to be tested. The chamber dimensions 
are equivalent to those of the ascent engine, l.e. , 7.79 inch diameter, con- 
traction ratio of 2.90 and L* » 27.6 inches. Stability rating was accomplished 
by disturbing the combustion process with 6.5 grain (RDX) bombs, with two 
bombs being used for each hot firing. 

HARDWARE DESCRIPTION 

The engine used for the hot firing program was patterned closely after the LM 
ascent engine with maximum use being made of hardware fabricated and used during 
the_preceding program. An assembly drawing of the engine is shown in Fig. 39; 
a photograph of the engine components is shown in Fig. kO. Provisions were 
made in the engine for installing instrumentation at a number of locations: 
eight bosses for Photocon high-frequency transducers, four bosses for Kistler 
high-frequency transducers, two bomb locations, one chamber pressure tap, and 
twelve thermocouple locations for cavity temperature measurements. The solid 
wall thrust chamber was constructed from the previously used chamber by removing 
the upstream portion of the chamber, including the pulse gun parts, and welding a 
new section in its place. The various components of the thrust chamber assembly 
and their function are summarized in Table 5» Further, the_thrust chamber 
instrumentation locations are listed in Table 6. 

The only hardware design and modification that caused any difficulty was the 
thrust chamber weld. This weld cracked after a number of tests. However, 
after re-weldlng, no further difficulties were encountered* 
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THRUST CHAMBER 



Figure 39. Full Scale Thrust Chamber 




Photograph of Full Scale Engine Components, 

Including Cavity Rings and Feed Line Quarterwave Tubes 
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TABLE 5 


PULL SCALE TEST HARDWARE 


ASSEMBLY 

SUB- ASSEMBLY 

DRAWING NO. 

COMMENTS 

Thrust 

Chamber 

LM Thrust Chamber 
Rebuild Assembly 

R3J1813 

Modified original to meet 
program requirements. 


Special Photocon Adapter 
Pitting 

R3B1814 

Allows a 1345-series photocon 
to be used in a bomb pox t. 


Kistler Transducer Adapter 

R3C1818-21 



6.5 Grain Bomb Assy 

R3C1823 

Drawing used for fabrication. 


Bomb Adapter Assy 

R3C1821 

For mounting the low cost bomb 
in the thrust chamber 

Adapter 

LM Baffles Injector Adapter 

R3E1815 

Adapts R3DI816A injector to 
R3J1813 thrust chamber. 


Ki stler Transducer Adapter 

R3Cl8l8-ll 



1 Kistler Transducer Adapter 
i Housing 

R3C1819 

Allows the cavity transducer 
position to be varied from 1.0 
to 1.8 inches deep. 


Adapter Housing Clamp 

R3C1820 

Secures R3C1819 Adapter to 
R3E1815 adapter. 


Injector Retention Plate 

R3C1817 

Secures the R3D1816A injector to 
the R3E1815 adapter. 

Injector 

LM Unbaffled Injector Assy 
Revork 

R3D1816A 

1 

Remove the flange from the 
original; this allows for radial 
cavity thermocouples to pass 
the. injector seal. 

Engine 

Assembly 


R3J-1822A 

Shows complete assembly and 
instrumentation location. 
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TABLE 6 


THRUST CHAMBER INSTRUMENTATION LOCATIONS 



TRANSDUCER 


High Frequency Chamber Pressure Photocon 


Time-Average Chamber Pressure Tabor 

High Frequency Cavity Pressure Kistler 


Radial-jCavity Temperature 


Axial Cavity Temperature. 


r 

(Inches ) 


3-895 


3.895 

3.895 

( 2 ) 

3.895 

( 2 ) 

3.395 

4.095 
4.195 
4.495 
* 1.795 

5.095 

3.995 


LOCATION^ 1 

0 

(Degrees) 

Z 

(Inches ) 

0 

4, 

.23 

270 

4, 

.23 

225 

4, 

.23 

135 

4, 

.23 

90 

4, 

.23 

0 

2, 

.06 

270 

2, 

.06 

225 

2, 

.06 

330 

5. 

.8 

330 

1 , 

.60 

330 

(2) 

30 

1 1, 

.60 

30 

(2) 

90 

(3) 

95 



85 



90 



95 



85 



120 

-0, 

.1 

125 

-0, 

.2 

115 

-0, 

• 3 

120 

-0, 

.6 

125 

-0, 

.9 

115 

-1, 

.2 


(1) The coordinate system is referenced to the center of the injector face, with 
the positive axial direction taken as downstream and the angular locations are 
measured from the vertical direction at the top of the engine. 

(2) r & 2 varied with the cavity depth and location. 

(3) thermocouple immersion varied with cavity loeation. 
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Thrust Chamber 


The thrust chamber used during the previous program was modified by replacing 

the upstream portion of the chamber vith a new section which was electron-beam 
welded in place. This weld cracked after about thirty hot firings. The crack 
began on the inner wall/ slowly propagated through to the outside and then 
began to leak. The external Joint was then successfully rewelded (tungsten/ 
inert gas). 

Special adapter fittings were designed for the bomb installation ports in the 
thrust chamber. Local deformation (splay) of the bomb ports had been encountered 
previously from the bomb detonation. The bomb adapter fitting was fabricated 
from 43UO steel and subsequently heat treated to produce a yield strength in 
excess of 200,000 psi.- These fittings have prevented deformation of the bomb 
ports throughout the test program, in excess of 100 bomb detonations per fitting. 

Injector Adapters 

Two injector adapters (R3E-I815) were fabricated; the first to accommodate 
both axial and radial cavity configurations,, and the second to accommodate 
only the radial cavity configurations. Two adapters were made so that the 
cavity thermocouples could be positioned before assembly. These adapters con- 
tain all of the cavity instrumentation ports (for thermocouples and Kistler 
transducers). In addition, they included the sealing surfaces for piston- type 
seals on the thrust chamber and injector. The. Kistler adapters were essentially 
an eccentric which allowed the transducer to be positioned at eavlty depths 
varying from 1.10 inches to 1.80 inches. „ 
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The injector adapters were designed- to eliminate numerous sealing problems 
encountered JLn the previous program. Pew problems were encountered with 
them. 

Injector 

The unbaffled injector used for the hot firing program was the same as that 
used in the previous program for the flight-design injector tests. This in- 

jector was made entirely from aluminum and consisted of a one-piece milled 
body, eight propellant rings and six propellant -distribution rings, which 
were electron-beam welded together. The propellant distribution rings were 
incorporated into the outer six manifold ring grooves to improve propellant 
distribution. 

Unlike -doublet injection elements were used in the injection pattern which 
was similar to the IM ascent engine injectors. All of the injection orifices 
were electrical-discharge machined. Standard IM ascent injector inlet screens 
were installed in the propellant inlets. 

The injection patterr was comprised of 21b unlike-doublet elements and. 30 axial- 
ly directed showerhead film coolant orifices. The impingement angles were 
28.5° on the oxidizer stream and k0° on the fuel stream. The. .nominal impinge- 
ment distance was 0.175 inches. Sixty of unlike doublets were comprised of 
0.0417-inch diameter fuel orifices and 0.0l40-inch diameter oxidizer orifices; 
these were used along the chamber wall. The remaining 15^ elements were com- 
prised of 0.0310-inch diameter fuel orifices and 0.0365-inch diameter oxidizer 
orifices. 


Ill 
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Two photographs of the injector, taken after completion of the hot -firing 
program, are shown in Fig. 41 and 42. The weld metal evident on. the per- 
iphery of the injector face- was added during the film cooling test series to 
close alternate orifices. The injector was in good condition in spite of 
the fact that it had been exposed to nearly 300 bomb disturbances. The in- 
\ jector was modified from its original design to remove a flange and substi- 

tute a piston-type 0-ring sealing arrangement. 

In the course, of the hot firing program only one minor injector repair was 
necessary. A small crack developed at the center of the injector face after 
a small protruberance- (the center disc) burned off. This crack was welded 
and no further problems were- experienced. 

Acoustic Cavity Rings 

Seventeen acoustic cavity rings were fabricated; these are summarized in 
Table 7 . The cavities were partitioned at 30° intervals. The partitions 
were -tng+.wneri ei ther wit h roll. pins or by welding. When installed, an 
adhesive (Permatex No. 2) was used to fill any unwanted cracks. A photograph 
of 15 of the cavity configurations ia. shown in Fig. 44 , subsequent to testing. 
rp^Q 0 f -^he original rings were modified to make the additional cavity 
configurations. 
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Figure Ul. Photograph of Injector After Completion 
of Test Program 



Figure h2. Photograph of Injector After Completion 
of Testing 


TABLE 7 
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(Cavity Dimensions Exaggerated 
for Clarity) 



Figure 1*3. Notation For Catity Configurations Evaluated 
During Phase I Hot Firings 



t 

< 

1 

I 

| 


116 





Stability Rating Bombs 

jj . 

!$’■ Bombs were used exclusively during the hot firing program to stability rate 

the engine. In addition, a new low cost comb configuration was developed and 
used throughout. A schematic diagram of the new bomb is shown in Fig, 45a j 
a photograph is shown in Fig, 45b. Although some development was. required, 
a substantial cost saving was achieved. The new bombs were obtained for 
less than $10 per bomb whereas the bombs used previously, throughout the IM 
ascent engine development and the acoustic cavity study, cost approximately 
$100 per bomb. Although the price for an additional purchase of the low-cost 
bombs would probably increase somewhat, the cost saving would still be 
substantial. 

The performance of the low cost bombs, in terms, of overpressure -and pressure- 
rise rate, appears entirely equivalent to that of the old design. Typical 
bomb disturbances obtained in the engine, one for an unstable test and one 
for a stable test, are shown in Fig 46 , Furthermore, the recorded over- 
pressures and pressure rise rates obtained throughout the hot firing program 
consistently exceeded the contract requirements, i.e., peak pressure in excess 
of 175 percent of the chamber pressure disturbance and e pressure rise rate 
in excess of 100,000 psi per second. 

Test Stand 

A schematic diagram is shown in Fig. 47 of the test stand and propellant sup- 
ply systems used for the hot firing program. Propellants utilized throughout 
the program were NgO^/^H^-UEMH (5O-5O). The run tanks were sufficient to 
permit cumulative run durations of 60 seconds. Gaseous nitrogen was used for 
propellant pressurization and system purges. 
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Run Ho. 57 
6.5 Grain Bomb 




With the Original Bomb 







Before- firings were begun, several modifications were made to the facilities 
to improve the control over engine mixture ratio and chamber pressure. Prin- 
cipally, this involved installation of cavitating venturies in both propellant 

lines and adding precision measurement and readout systems for each tank 

pressure. A digital voltmeter was used for the readout, the pressure signal 
being conditioned so that the readout could be interpreted directly as a pres- 
sure in psi. 


Measurements 


The instrumentation has been summarized previously in Table 6 . Conventional 
measurements of pressures, propellant temperatures and flowrates were made to 
determine the steady state operation. Cavity temperature measurements were 
made with tungsten/rhenium thermocouple (W26R e /W5R e ); these thermocouples will 
withstand temperatures in excess of 3000 F. High frequency pressure measure- 
ments were made with Photocon transducers (model I 3 U 5 ) and Kistler transducers 
(model 603A), the latter being using to measure the cavity response. The 
steady state operational data and the temperature data were recorded on a Beck- 
man 210 Data Acquisition System. This is an analog-to-digital conversion sys- 
tem with an established accuracy of <0.1 percent. The high frequency pressure 
data were recorded on a 14 channel Hewlett-Packard model 3955 ® tape recorder. 
No significant problems were encountered with any of the measurements during 


the program. 


PRASE I - FULL SCALE TEST SERIES 

Test Plan _ 

A series of 45 hot firings was planned for the Phase I full-scale test program; 
however, a number of additional tests were required because of operational and 
hardware difficulties experienced during this series. These tests were directed 
toward experimental definition of the effects of cavity width, axial location 
and multiple rows on stability. In addition, the cavity temperature distri- 
bution and oscillatory pressures were measured. 

Several sets of tests were planned for this series* At the_outset several 
tests were planned for reference purposes and to gain assurance that the sta- 
bility characteristics of the engine have not changed significantly from those 
observed during the IM cavity program (Ref. l). A 1.3-inch long, axially dir- 
ected cavity was selected; this cavity length had previously produced marginal 
stability when the width was 0.2-inches and complete stability when 0.3-inches 
wide. These tests were expected to be sufficient to determine any changes in 
stability. 

The next set of tests was made to determine the effect of cavity width on sta- 
bility. During previous testing a. minimum effective cavity width was found 
(Ref. l), but no maximum width. Therefore, a set of tests was -made to deter- 
mine this upper limit. Cavity widths of 0.3, 0.6, and 1.0-inehes were 
evaluated. These tests Were expected to be adequate to define the upper 
bound on width. 
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The third set of tests was planned to investigate the effect of axial position 
on cavity stabilization. Previous tests had indicated a strong sensitivity 
of cavity effectiveness on axial position (Ref. l). These tests were planned 
to verify this effect and/or examine it further. Specifically, tests were 
scheduled to Investigate the variation of cavity effectiveness with axial 
position and, further, to determine if the previously observed Iobs in effect- 
iveness is only due to a change in required cavity tuning. The intent was to 
determine the criticality of cavity position. 

The final set of tests was planned to investigate the potential advantages of 
multiple rows of cavities. The intent of these tests was to provide a measure 
of the ability to increase cavity damping by adding resonators. 

All of the Phase I hot firing were to be made at nominal IM ascent engine 
operating conditions, i.e., a mixture ratio of 1.60 and a chamber pressure 
of 123 psia. Propellants were (green) and (50-50) at ambient 

temperature. Facility modifications necessary to restrict mixture ratio 
variations to within 1.60 ± 0.03 and chamber pressure variations to within 
123 ± 3 psi, on a 2 a basis, were made. 

Bombs (6.5 grain) were used exclusively for stability rating. 

Operational and Stability Results 

A total of fifteen cavity configurations was tested during this test series* 
These configurations, the corresponding test objectives, and the observed 
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engine stability are summarized in Table 8 , the completed test matrix for 
this test series. The engine was considered stable if it damped at least 
six bomb disturbances, generally with two being used for each hot firing. 
Clearly the engine was rather effectively stabilized by most of the config- 
urations tested. The operational and stability results, along with the char 
acteristics of the bomb disturbances, are summarized in Table 9 . The bomb 
over-pressures and rise rates fall within an acceptable range j they are not 
substantially different from those obtained previously under the IM ascent 
engine contract with the more expensive bomb configuration. Further, the 
variations in mixture ratio and chamber pressure are also within acceptable 
limits. The mean values of chamber pressure and mixture ratio are 122.47 
psia and 1..608, respectively. Arithmetic average deviations on chamber pres 
sure and mixture ratio arejQ.819 psia and 0.010 5, respectively. Root -mean- 
square deviations (standard deviations) are 0.967 and 0.013. The cavity 
temperature data have been found to be repeatable and quite consistent. 

These data will be discussed later in this section of the report. 

Several important observations may be drawn from the stability results . 
First, no slot width was found which is "too" -wide; the maximum effective 
slot width is greater than 1.0 inches. Second, the engine can be made sta- 
ble with a slot located as far as 1.0 inches (midwidth) downstream if the 
slot is wide enough. Finally, with the dual axial/radial configuration (No, 
14), although the engine was stable, the damp times were somewhat longer 
than with other configurations, suggesting a smaller stability margin. 
However, this probably reflects the same error in hardware design which 
caused the axial slot results to be deemed invalid , as discussed below. 










1.35 axial .2 1.3 0 122.4 1.615 1.1 325OHZ OSC. 
1.2 I I I I 121.2 1.59 .9 325OHZ OSC. 
1.2 i T f r 122.8 1.609 .9 3400HZ OSC. 




SUMMARY OF PHASE I FULL-SCALE TEST RESULTS 
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The stability results obtained with the 0.2 and 0.3-inch wide axial slots 
were invalid because the slot entrance was partially blocked by the radial 
blank ring. The blank ring which fills the radial cavity location must be 
chamfered to completely expose the axial cavity; it was not. As assembled 
the exposed slot width was only ~ 0.1 inch, which is too small to stabilize .. 
the engine even if it were not detuned by the obstruction. 

A number of hardware and operational difficulties were encountered which — 
prevented the efficient completion of this test series. Seme of the oper- 
ational difficulties stemmed from the facility modification made to improve 
control over the propellant flowrates and pressures. One problem, that may 
have affected the stability results in some small measure was a persistent, 
low-level feed-system coupled oscillation at 500 Hz that was repeatedly ex- 
perienced during the testing. The worst levels were ~ 10 percent of chamber 
pressure but generally it occurred at ~ 5 percent peak to peak. Initially, 
the cavitating venturies were moved well upstream of the injector ( > 160 
diameters) to eliminate the possibility of two-phase flow in the lines. When 
the problem persisted, quarterwave tubes were added to each feed line a few 
inches upstream of the manifold in an effort- to reduce or eliminate the 
problem. Each of these consisted of a 1-inch-diameter tube mounted vertical- 
ly above the feed line and coupled through art orifice to another closed length 
The upper cavity is intended to accept all of the entrapped gas as the line 
primes. The liquid-filled tube to the orifice should then act as a quarter- 
wave tube. This approach allows the quarterwave tubes to drain naturally 
without the need for bleed valves. These tubes were installed subsequent to 
test No, 6k . Tests 65 through 7 6 were made to repeat previous tests with 
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a 500 Hz oscillation level ~ 10 percent. The levels were reduced by the 
quarterwave tubes to ~ 6 percent. 

The influence of the 500 Hz oscillation on the high frequency stability of the 
engine is difficult to adequately assess. However, there is no evident difference 
in the bomb responses observed when the amplitude of the 500 Hz oscillation 
varied between 5 and 10 percent. Moreover, this amplitude range is low enough 
that mode coupling, between this oscillation and any high frequency mode, is 
not ejected to be important. Thus, these levels are believed low enough so 
that stability of the acoustic modes is substantially unaffected by the feed 
system-coupled oscillation. However, the cavity temperature may be increased 
somewhat so that the required cavity tuning may be modified. 

A similar 500 Hz oscillation was experienced during the LM ascent engine injector 
development, as reported in phase reports for both the Design Feasibility Itest 
Ehase and the Design Verification Test Phase. With the unlike-doublet injector 
the problem was eliminated by changes to the fuel duct, the changes being 
intended to reduce start transient over pressures. Use of this duct during 
previous testing with the unbaffled injector under the LM engine contract may 
be the reason_lhat no $00 Hz oscillation was experienced during that testing. 


OFF- NOMINAL FULL-SCALE TEST SERIES 
Test Plan 

A series of 15 hot firings was planned for the off-nominal test series to 
evaluate the influence of engine operating conditions on cavity effectiveness 
(stabilization). A single cavity configuration was selected which was ex- 
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pected to be near the stability limit, but adequate to stabilize the engine 

at nominal conditions; this configuration was a 0. 3- Inch wide, 1.3-lnch deep, 

axially directed cavity. A near-marginal configuration was selected so that 
changes In the stability of the engine would be more nearly evident. 

The off-nominal testing was Intended to show the sensitivity of cavity 
stabilization to changes in the operating conditions. A strong sensitivity 
was not expected, based on results from the 1M ascent engine contract. How- 
ever, it was deemed important to define the sensitivity of cavity stabiliza- 
tion to operating conditions because, in particular engine applications, close 
control of engine operating conditions may not be possible or even desirable. 
Thus, the applicability of cavities could be determined in part by this degree 
of sensitivity .to operating conditions. 

One reason for selecting the cavity configuration indicated above was to 
re-evaluate it because of the failure to obtain valid re stilts during the 
Phase I testing. During the Phase I tests this cavity configuration failed 
to stabilize the engine because the cavity entrance was partially blocked by 
the improperly machined cavity rings. However, this configuration proved 
adequate during the LM engine program (Ref. l) with the same injector, there- 
fore, it was expected to be adequate and appropriate for the planned tests. 

A test matrix for the planned tests is shown in Table 10. One set of tests 
was planned at nominal conditions to verify stability. Subsequently,, the 
effects of first chamber pressure and then mixture ratio- were to be determined 
sequentially by variations from the nominal. A minimum pressure of 90 psia 
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TEST MATRIX FOR OFF-NOMINAL FULL-SCALE TEST SERIES 



was selected as probably the lowest allowable without introducing significant 
chugging instability. A maximum chamber pressure of 200 psia.vas selected. This 
pressure should be high enough to be significant but would not require a facility 
change. The mixture ratio range was selected arbitrarily but was believed large 
enough to produce a significant effect on stability, if stability was affected 
by it. 

Operational and Stability Results 

Stability results obtained from this series of tests are summarized in Table 11, 
the completed test matrix. The operational results and bomb response character- 
istics are summarized in Table 12. During this series, fourteen tests were 
made with a 0.3-inch by 1. 3-inch axial cavity with chamber pressure and mixture 
ratio being varied from nominal conditions. Only two firings, rather than 
three, were made at the high mixture-ratio condition because scheduling problems 
in the test area made the third test expensive. Actual test conditions were 
slightly different from planned; the tests were made at chamber pressures of 
190, 120 and 90 psia and mixture ratios of 1.41, 1.60 and I.78. 

It should be noted that the tests at nominal conditions were unstable during 
this test series, whereas previous tests with the same cavity configuration 
were stable during testing under the IA ascent engine contract (Ref. l). Oils 
difference may be attributed to moderate chugging (the amplitude is approximately 
7 percent peak-to-peak of chamber pressure at a frequency of 500 Hz) observed 
during the current test series. The chugging mode could alter the gas tempera- 
ture in the cavity which, in turn, would alter the cavity tuning and damping. 


Y OF OFF-NOMINAL TEST RESULTS 
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The test results show that there is a definite effect of operational conditions 
on engine stability. A gradual worsening of stability is indicated with in- 
creasing chamber pressure or decreasing mixture ratio. However, these results 
exhibit only the overall effect on stability; it is not poss ible to separate 
effects on cavity damping from those on the net driving processes. Nonetheless, 
clearly there is no abrupt or drastic change with operating conditions. The 
observed effects are probably due to changes in the combustion driving rather 
than changes in cavity damping. It does not appear that cavity damping is 
seriously degraded by changes in operating conditions. 

FILM-COOLING FULL-SCALE TEST SERIES 
Test Plan 

Hie effects of varying the film-coolant flow- rate on. cavity stabilization were 
evaluated in a series of tests because detrimental effects of this cooling on 
stability have been observed with acoustic liners (Ref. 13). One possible explan- 
ation for this observation is that part of the film coolant is trapped in the 
cavities which would, alter the tuning and damping of the cavities. Another possi- 
bility is that the oxidizer may penetrate into the film-coolant boundary layer and 
undergo a chemical reaction. In the latter case, the liner surface could act 
effectively as an acoustic energy source rather than an energy sink. Therefore, it 
was anticipated that film cooling could also pay an important role in stabiliza- 
tion with acoustic cavities. 

A series of 18 hot firings was planned to investigate these effects^— _Filra- 
coolant flow-rates of 50, 100 and 200 percent of nominal were, planned, where 
the film coolant flowrate used on previous tests was regarded -as 50 percent 

rather than 100 percent of nominal. — This choice allowed the flowrate to be 
varied by drilling additional equally spaced coolant oriflces-and then 
selectively plugging them. Further, the 100 percent flowrate is near the 
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nominal film coolant flowrate on the production IM ascent engine injector. 

The selected nominal flow-rate is roughly 4.6 percent of the fuel flow-rate 
(this percentage is based on total orifice area rather than measured 
flow-rates) . The other operating conditions were to be kept at their 
nominal values. 

Two cavity configurations were planned this series: (1) a 0.3-inch wide 

1.3-inch deep, axially -directed cavity, and (2) a 0.4-inch wide, 1.5-inch 
deep radially -directed cavity. The 0^3 x 1.3-lnch axially -directed cavity 
had been shown to be a marginally unstable configuration (off -nominal., test 
results) and the 0.4 x 1.5 -inch. radially-directed cavity was a stable con- 
figuration (Phase I full-scale test results). These two cavity configurations 
were selected to provide a stability regime such that the sensitivity of 
acoustic cavities to film-coolant rate can easily be distinguished. 

The injector was modified to achieve the experimental objectives. Initially 
the injector had thirty #80 eloxed film coolant orifices. To these, 90 addi- 
tional holes-were. added to provide 120 equally spaced #80 holes required for 
the 200 percent film-coolant flowrate. The 100 percent fllm-coc\ant rate was 
to be achieved by closing alternate holes and the 50# level was then to be achieved 
by closing every other remaining hole. The orifices wore welded closed. 


Operational and Stability Results 

A test matrix for the planned and completed tests is shown in Table 13* 

The observed engine stability is also indicated there. The operational 
results and observed bomb response characteristics are summarized in Table lit. 
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OF FIIM COOLING FULL SCAIE TEST RESULTS 
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Nineteen full-scale tests were made for the film-cooling test series. Two 
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cavity configurations were tested at each film-coolant flow-rate# Mixture 
ratio and chamber pressure for these tests were adequately controlled, mixture 
ratio being maintained at 1.60 + 0.03 (o/f) and chamber pressure at 126 ± 2 psia. 


Examination of these. test results clearly shows that the film coolant flowrate 
has little, if any, influence on stabilization with acoustic cavities. No 
doubt, film coolant characteristics can influence the cavity temperatures 
(sound velocity) , and thus tuning, under some circumstances. But, the stabi- 
lizing ability of the cavity is not degraded, as opposed to the. reported effects 
with acoustic liners. Therefore, fortunately acoustic cavities may. be designed 
without careful tailoring of the film coolant characteristics simultaneously. 
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UNCONVENTIONAL- CAVITY FULL-SCALE TEST SERIES 

The test plan and cavity design for this test series was approached somewhat 
differently from the other full-scale test series. In part, this was done 
because a series of tests was originally planned to verify the effectiveness 
of the developed design methods. Subsequently, this demonstration was deleted 
in favor of the unconventional cavity series because the latter appeared more 
valuable. However, because the original purpose of the test series was one of 
verification, the cavity selection for the unconventional test series was made in 
a manner similar, insofar as possible, to that which would be used for a demon- 
stration of the design methods. Thus, the intent was largely to design and 
demonstrate stable absorber configurations, rather than to survey effects by 
testing a range of configurations. 


3M 




A series of 12 hot firings was planned to evaluate the damping capability and 
limitations of unconventional acoustic cavities. Such cavities are simply 
acoustic resonator configurations which differ from those in common use for 
stabilization, i.e., Helmholtz or-quarterwave resonators. Two kinds of 
unconventional cavities were to be tested; these kinds have been denoted 
as "intermediate" and "L-shaped" resonators. Analytical and subscale experi- 
mental studies done earlier, in the_prpgram were used to aid in characterizing 
these resonators. 

The principal .reason for interest in unconventional cavity configurations is 
to add versatility and latitude in cavity design so that maximum use of the 
available space in an engine can be made. Furthermore, the analytical studies 
Indicate that some damping advantages can also be realized, e.g., greater 
bandwidth and multiple resonant frequencies. 

A test matrix for the planned tests is shown in Table 15 . Four cavity con- 
figurations were selected for evaluation and these are shown in the matrix. 

Three hot firings, with two bomb disturbances (6.5 grain bombs) introduced 
during each hot firing, were planned with each cavity configuration. The 
tests were to be made at nominal operating conditions, i.e., a mixture ratio 
of I.60 and a chamber pressure of 122 peia. 

The hardware designs for the four cavity configurations to be tested are 
shown in Fig. k8 (except the Kistler transducers were omitted) . All cavities 


TABLE 15 

PLANNED TEST MATRIX FOR UNCONVENTIONAL CAVITY 
FULL-SCALE TEST SERIES 


Test 

Objectives 

Aperture 


Cavity 

Number 

of 

Width (in.) 

Length (in.) 

Width (in.) 

Length (in.) 

Firings 
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had radially directed entrances and were machined into a ring which fit into 
the existing thrust chamber, without modification. That is, the cavities 
were designed to fit within the existing 1.3 x 1.8-inch envelope. Suitable 
cavity configurations with axially directed entrances could not be made 
without modifications to the thrust chamber and, therefore, were not 
selected. 

These cavity designs were obtained from damping calculations based on the 
subscale (T-burner) test data and the previous hot firing results. One 
"L" shaped and one intermediate cavity (the larger of each kind) were 
designed for an average cavity temperature of 1900°F, which is the average 
temperature measured during the Phase I full-scale tests with 0.3- to 0.4- 
inch-wide radial slots.. However, this temperature is only an integrated 
(linear) average from the thermocouple data; therefore, it may not be an 

appropriate average or effective value for the acoustic calculations. 

Consequently, a second pair of cavities vs s designed by referencing the 
Phase I full-scale stability results, which suggest a maximum damping (or 
at least adequate damping) at a slot length of 1.5 inch. An effective 
sound velocity corresponding to maximum damping with a quarterwave slot 
at 1.5 inch was, therefore, used to design the second pair of cavities. 

This effective sound velocity corresponds to an average cavity temperature 
of 900 a F . Predicted damping curves for each of these cases is shown in 
Fig. 49 through 52 . These curves were obtained with the area-weighted 
cavity admittance (before the error in the computer program was-discovered) . 
The cavity dimensions were chosen for maximum damping. 
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Figure 50. Predicted Damping for "L" -Shaped Cavity at 900F 
(With Area-Weighted Cavity Admittance) 
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The results from the T-burner tests enter the design calculations through the. 
resistance coefficient, r . However, the experimental results yielded values 
of r which appeared higher than reasonable. Consequently, somewhat lower 
values were used for the design calculations; values of 1.5 and 2.0 were used 
for the intermediate and "L" shaped cases, respectively. The corresponding 
values obtained from the T-bumer data were. 1.8 and 4. Although the amount 
of this reduction is rather arbitrary, other damping calculations_indicate 
the effect of changing T in this case is principally to vary the magnitude 
of the predicted damping but not the cavity dimensions corresponding to 
maximum damping (at lea st fo r practical purposes). Thus, the same designs 
would probably be selected even if the high values of T were used. 


Operational and Stability Results 

The completed test matrix for this series of tests is shown in Table 16; all 
tests were stable. The operational results and bomb response characteristics 
are summarized in Table 17. Thirteen hot firings were made during which 
the mixture ratio was 1.60;fc 0.03 and chamber pressure was 123 i 2 psia. 

These results clearly demonstrate that unconventional cavities can be effec- 
tively used in place of simple Helmholtz or quarterwave resonators. Further- 
more, adequate stability was obtained with a relatively large variation in 
cavity tuning, suggesting either a high margin of stability or broad bandwidth 
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CAVITY TEMPERATURE MEASUREMENTS 


Relatively detailed temperature measurements were made during essentially 
all hot firings. Tungsten rhenium thermocouples (W26 Re/ws ft® ) were used 
for these measurements because of the ability to withstand relatively high 
temperatures. During the LM ascent engine acoustic cavity study (Ref. l), 
cavity temperature measurements were made but these were not of sufficient 
detail to accurately define the temperature distribution within a single 
cavity. Therefore, the measurements during the current program were directed 
toward defining this distribution. 

Exposed Junction thermocouples were located in a radial cavity as shown in 
Fig. 53* Six thermocouples were used, these being 0.1, 0.2, 0.3, 0.6, 0.9 
and 1.2 inches from the open end of the cavity. The three thermocouples 
nearest the open end were made with 22 gauge (0.025 inch diameter) wire, 
butt welded to form the junction. The remaining three were made of 30- 
gauge (0.010-inch-diameter) wire, also butt welded to form the junction. 

With a radial cavity , the thermocouples are located around an angular posi- 
tion of 90 degrees from the vertical (top), the thermocouples being staggered 
+5 degrees to give adequate spacing as shown in Fig. 53. With an axial 
cavity, a similar arrangement was used but the thermocouples are located 
around 120 degrees. The fuel inlet was located at 270 degrees. 

The temperature data obtained in this manner were found to be repeatable 
and quite consistent. Extensive data were obtained; these are tabulated in 
Table B-l of Appendix B* Some representative examples are shown in Fig. 54 
and Fig. 55. In addition, art attempt was made to develop or adapt an 
analytical model to describe the temperature distribution. 
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Figure 5^. typical Cavity Temperature Data from 
Phase I Full Scale Tests 
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This analysis is described in Appendix B. Unfortunately, the effort was not 
entirely successful; however, a semi -empirical equation was developed from 
the data which may be used for rational estimates of the cavity temperature 
distributions. This analysis suggests that cavity temperatures are most 
strongly related to: (1) film coolant characteristics, (2) local adiabatic 
wall temperatures, and (3) near-wall propellant injection characteristics. 

Naturally, the propellant combination Itself affectB the temperature distri- 
bution as well. 

Average cavity temperatures, obtained by numerically integrating the data along the 
slot length, are shown in Fig. 56 for the Phase I tests. The effects of axial 
position and slot width on this average temperature are clearly shown. Similar aver- 
age temperatures are shown in Tables 18 through 20 corresponding to the data 
obtained during the other test series* The latter results are somewhat 
less consistent but some trends are evident. The average cavity temperature 
appears to increase with chamber pressure and possibly show a maximum, but 
the latter trend is dependent upon the low indicated temperature at high 
mixture ratio which may be suspect* The average temperature results in 
Table 19 suggest an increase with decreasing film coolant flow rate, as 
might be expected. The results from the unconventional cavities. Table 20 , 
appear consistent with similar conventional cavities, suggesting that the 
cavity shape does not have a substantial effect* The thermocouples for 
this test series were located near the midwidth position for each radial 
location* 
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Spatial Average Cavity Temperature Measured 
During Phase I Full Scale Tests 
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table 18 

iVUnaoE CAVITY TEMPERATURES MEASURED EURIEO 
SPATIAL AVERSE CAHT^^ ^ 


Mixture Ratio 

(0/F), 

1.60 

1.60 

1.60 

1.41 

1.78 


Average Cavity- 

Temperature (F]_ 

1440 

1200 

1750 

1330 

780 


table 19 

SPATIAL AVERAGE CAVITY lll * IM0 

SPAii FIIM COOLING test series 


Coolant Flowrate 
Percent of Nominal 


A ^^Temperature (F) 

0.4 x 1.60 in. Radial 0-5 * *** 

Cavity - — 


1530 

1440 

1430 


1770 

1370 

1540 


TABLE 20 

■“JSSJ5S.” ”” 


Configuration 


Intermediate 
Intermediate 
"L" -Shaped 
"L" -Shaped 


Resonator Cavity 
Width (in.) 


Average Cavity 
Te mperature (Fl 

1940 

1900 

1980 

1780 



The foregoing results indicate that the cavity temperature depends in a 
rather complicated way on cavity width and axial position and on engine 
operating conditions and film coolant characteristics. Moreover, a_ com- 
plicated dependence on slot depth is also indicated. Hiese results can be 
introduced into the damping calculation to the extent that the complicated 
temperature profiles can be approximated by a smooth curve .or a uniform 
average temperature. A more sophisticated analysis of the cavity wave motion 
which accounts for the details of the temperature distribution is probably 
unwarranted at this stage because of uncertainties in other factors such as 

molecular weight of the gas. Nonetheless, the average values are required 

for cavity design and are probably sufficient. 


SUMMARY OF FULL SCALE TEST RESULTS ...... 

Results from the full scale testing allows clarification of several aspects 
of cavity behavior. The detailed temperature data just described appear more 
precise and consistant than any previously available. These data can be 
averaged with greater confidence than previously. Moreover, the influence 
on cavity temperature of cavity width and depth, and of engine operating 
conditions and film coolant flowrate are reasonably clear from the data. 

Further, the data from the unconventional cavities suggest the temperatures 
therein are largely equivalent to those in conventional cavities. 

\ 

During this program a low-cost bomb was developed which exhibits the same 
yield as previous bombs but with a substantial cost saving. 

\ 
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Puri* the Phase I «U »««*. th. Wldth ’““ 

location and multiplicity on engine etabiUty were meaeured. B» re. 

indicate a relatively Oarge range of oavity di^neion. ewiet. which w 

effect .lability. Ho clot width up to 1.0 inch., wae found "too w * ° 

stabilize the engine, the maxima effective .lot width wa. not reachM-ing 

.. anpine was readily stabilised with radially 
the program. In addition, the engine was 

directed .lot. local* up to 1.0 inches (nidw^th) downstream *- 
face, thus, th. stabilization ability of the slot is not highly sensit v. to 
position a. some previous data had suggested. Moreover, the engine wa. 
siso stabilized with two and three .lot combination., however, the o^ 

evident advantage of -*U — - * — " ^ ^ ^er- 
srsa is to introduce multiple resonant fmuencie. of the cavities. Furl 

M r„ during these test. the influx, of engine operating options was 

^ by precis, control of the engine mixture ratio - chamber pr*.ure. 

The influence of operating conditions on cavity stabilisation was 

4 a The test results exhibit a gradual and 
during th. off ncrfsml t«t »” iM - tM 

— — •• - " :. irxr™ 

processes rather than a cavity effect, me s .a 

does not appear highly sensitive to engine operating condition.. 

„ Influence of the film coolant flow rate on oavity stabilization was a.o 

investigated, m. etability results indicate little influence, if • 

.v. range investigated, me oavity tea^eratures 
film coolant flowrate over th. range inveaerg 

J . hot a slisht change of etability with the 0.3- 
were affected to some degree but a eiignz cmuik 

inch wide cavity appears edible With th. temperature change, m. observe 
wash dependence of cavity stabilization on film coolant flowrate remove, con 
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cem for a degradation of stability accompanying increased film coolant 
flowra te which has been observed with acoustic liners. 

The stabilization ability of four different unconventional cavity configura- 
tions (two "L" shaped and two intermediate configurations) was evaluated. 

All four produced dynamic stability. This demonstrated effectiveness allows 
greater design flexibility to be introduced and more efficient utilization 
of available space. Moreover, some stability gains were achieved; the 0.3- 
inch-wide intermediate cavity stabilized the engine whereas none of the 0.3- 
inch-wide conventional cavities completely stabilized the engine. Further- 
more, the unconventional cavities were designed for two widely separated 
cavity temperatures and, thereby, resonant frequencies. Nevertheless, all 
of them stabilised the engine, a result which suggests a substantial stability 
margin. 

During the full scale testing the 0.3-inch wide, 1,30-inch deep axially 
directed cavity failed to completely stabilize the engine whereas during the 
114 ascent engine study the same configuration stabilized essentially the 
same engine. This difference may have been due to the use of different feed 
systems, which may have affected propellant spray atomisation and combustion 
characteristics (and thus stability), or due to the 500 HZ chug (or buzz), 
which could have altered cavity temperatures or high frequency stability. 

Perhaps another possibility is the difference in bombs used. Nonetheless, 
the stability differences are believed relatively small in that this configura- 
tion appears near the stability boundary from both, sets of stability data. It 
is reasonable to expect that near a stability boundary s m a l l changes in a 
variety of parameters could alter stability from being stable to being marginally 
stable as observed. 
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CURRENT DESIGN STATUS AND 


RECOMMENDED PROCEDURES 

The stability results from the Phase I Pull Scale Test Series may be con- 
veniently summarized in graphical form as shown in Fig, 57 . These results 
were- all obtained at nominal operating conditions. Clearly, a relatively 
large region of stable cavity configurations, in terms of cavity width and 
depth, exists for this engine. However, the boundaries of that region have 
not yet been defined, except for the smaller cavity width. The indicated 
size of this stable region suggests the ability to stabilize more unstable 
engines with acoustic cavities than the engine tested, which was, nontheless, 
quite unstable without a cavity. 

Also shown on this plot are predicted curves of constant damping coefficient. 
These curves were calculated for radially directed slots located adjacent to 
the injector and with an average cavity temperature of 900 P. The curves 
shown in Pig. 57 were obtained with rp/7P 0 - 0.18. Similar curves are 
shown in Fig, 58 for rp/7P 0 = 0.09. Of course, not all of the data were 
obtained with radial slots located at the injector as assumed in the calcula- 
tion; nonetheless, the calculated behavior is generally in agreement with- 

the data in both cases. An accurate comparison between measured and calcu- 

lated damping is largely prevented because the boundaries of the stable 
region are not defined. In addition the measured average cavity temperatures 
varied with axial position and slot width so that a comparison with pre- 
dicted damping curves for a constant temperature is not quite valid. Further, 
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COMPARISON OF PREDICTED TRENDS WITH OBSERVED 
FULL SCALE STABILITY 
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the measured temperatures are significantly higher than that assumed for the 
damping calculation (£11900 F compared with the 900 F assumed in the calcula- 
tion). This difference may be due, in part, to the molecular weight of the 
gas in the cavity being higher than assumed (the main stream molecular weight) . 

In spite of these qualifying factors, the analytical model has several 
weaknesses. Although much effort has gone into assuring accurate prediction 
of the damping provided by an acoustic cavity, considerable uncertainty still 
exists. Furthermore, the model accounts for nonlinearities in only an approx- 
imate way, which must be regarded as a weakness. The importance of steady 
flow has not been demonstrated, but these effects can and should be included 
in the analysis. The temperature and sound velocity of the gases in the 
cavity cannot be accurately predicted. The damping model has been formulated 
to avoid analysis of the complicated combustion and other processes j however, 
this omission is also a weakness. 

On the other hand, the analytical model has been strengthened in a number, of 
ways. The wave equation is accurately solved for a chamber containing an 
acoustic cavity. The semiempirical cavity impedance expression that has 
been formulated and adopted appears largely satisfactory. Certainly, the 
ability to analyze the L'LH-shaped and generalized resonator configurations 
is important. 

From an experimental standpoint, considerable information has. been developed 
which is useful and important without regard to the damping analysis. A 
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wide range of cavity configurations has been found that will stabilize the 
engine at nominal and probably off-nominal operating conditions. These 
include simple quart erwave resonators, either axiall y or radially directed 
and located at the injector face or moved downstream. They also include 
the unconventional configurations. In addition, an extensive set of rela- 
tively accurate cavity temperature data has been obtained. 

RECOMMENDED DESIGN APPROACH 

Generally the recommended approach involves use of the damping-coefficient 
calculation to select a near optimum configuration within the available 
space which must contain the acoustic oavity. The optimum configuration 
may correspond to maximum predicted damping or seme compromise between high 
damping and weak sensitivity to poorly known or controllable parameters, 
such as cavity temperature or cavity dimensions. The iterative fora of the 
characteristic equation is recommended (either Eq. 35 or i±S ) without gains 
or losses (zero admittances) at the ends of the chamber. The generalized 
cavity impedance expression should be used whenever possible so that maximum 
use can be made of available volume. No preference is given for axially or 
radially directed cavities; either may be used. However, often greater open 
aroa within acceptable spatial limits can be obtained with a radially directed 
cavity. 

Selection of the cavity sound velocity and density should be based on what 
ever data are available for the propellant combination and injector type of 
interest. These parameters are generally difficult to select with confidence. 
For propellants, operating conditions, and injector types similar to those 
of the LM ascent engine, the cavity temperature data from this program may 
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be used as a guideline. Some of the available references for data on other 
propellant combinations are included in Appendix A. The cavity sound 
velocity and density may be controlled to sea* extent by purging the cavity 
with cold gas, e.g., hydrogen when it is also used as the fuel. Unfortunately 
the required purge flowrate is unknown. 

Based on the assumed cavity conditions, parametric curves of predicted 
imping as a function of, for example, cavity width or depth may be calculated 
from the characteristic equation. Gener all y, the zero-iteration-approxima- 
tion is sufficient, but this may be checked by adding iterations. Contour 
plots similar to those shown in Fig. 57 and 58 can be obtained by cross plot- 
ting. From these damping curves a suitable cavity configuration may be 
selected. Recognizing that many uncertainties are involved, it is not 
necessary to precisely define the point of optimum damping. S imil a r damp- 
ing curves may be calculated to exhibit the influence of cavity temperature; 
these affect the cavity selection to the extent the cavity temperature 
(sound velocity and density) are demaed uncertain. 

In making these calculations, several guidelines may be used. First, it is 
recommended, based on the most recent results from this program, that the 
slot amplitude/impedance parameter ffc/S P 0 he treated ae a combined parameter 
and that values for this parameter of 0.1 to 0*2 be used. This corresponds 
to selection of a design amplitude. Secondly, for inetability frequencies 
lower than say 1000 to 2000 Hz the volume requirements for a cavity become 
excessive and use of baffles rather than absorber* should be considered. 


The recommended chamber length for use in the calculation is the sum of the 
cylindrical length of the chamber and two-thirds of the distance from the 
beginning of nozzle convergence to the throat. Absorbers corresponding to 
damping coefficients less than ~ 100 sec”'*' are probably not worth hot fire 
evaluation. Partitions are recommended in the acoustic cavity to prevent the 
circumferential flow of hot gases. These should be spaced more closely than, 
perhaps, one-quarter wavelength based on the main chamber sound velocity. 

The cavities should be located near the periphery of the injector but the 
full-scale test data indicate that it is not critical that the cavity be 
immediately adjacent to the injector face. If more than one instability 
mode is of concern, then multiply tuned cavities may be used or the multiple 
resonant frequencies of the generalized resonators may be used. For example,., 
the lowest two resonant frequencies in a "L"-shaped resonator can be adjusted 
scL_that each corresponds to an instability mode. The full-s.cale_t.est data _ 
indicate that it is difficult to use a slot that is too wide; a result which 
suggests that preference be given to high open area designs. Further, these 
data indicate that film-coolant flowrate and engine operating conditions do 
not affect the cavity design except through the cavity sound velocity and 
temperature. 

Acoustic cavities may be constructed from ablative material as well as metals. 
However, in any case where dimensional stability is difficult to control, a 
design should be sought that is weakly sensitive to dimensional changes. The 
intermediate and generalized resonators are better in this respect that a 
straight slot. 

By employing this general approach and these guidelines it should be possible 
to repidly achieve a workable and stable configuration. 
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appendix a 


REVIEW OF PREVIOUS WORK 


INTRODUCTION 

Acoustic cavities or quarterwave resonators are part of a class of combustion 
stabilization devices which may be described as acoustic absorbers. The 
term "absorbers" is taken to delude all kinds of acoustic resonators with 
any distribution within the combustion chamber. The effects on stability 
of each kind of resonator array appears to be much the same in each case. 
Consequently, a review of the available information concerning the use of 
absorbing devices for combustion stabilization is the subject of this appendix. 

Various kinds of acoustic absorbers have been evaluated in a variety of 
engines. Much of the early work was directed toward the use of so-called 
full length acoustic liners. Such liners are comprised of an array of Helmholtz 
resonators distributed uniformly over the wall of the thrust chamber. In 
most cases the resonator array has been formed by installing a perforated 
shell with the perforations leading from the combustion chamber into an 
annular resonator cavity between the shell and the outer thrust chamber wall. 
The ability of full length liners to prevent high frequency combustion in- 
stability has been demonstrated. 

More recent work has shown that a full length liner is often not required 
for stability. Partial length liners were tried and found to be as effective 
as their full length counterparts. Still more recent work has shown that 
stability can be obtained with only a few acoustic resonators (often not of 
the Helmholtz type) located in or near the injector face* 
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Ahaorber fl nnflgurat*™" ayA location 

In general, acoustic absorbers are comprised of an array of acoustic resona- 
tors, i.c. Helwholtz, quarterwave, or other kinds of resonators. A Helmholtz 
resonator consist, of a small passage (an aperture) leading frm the main 
h amhS 2 lJLnto_a closed cavity (a resonator cavity). Early liner, were made up 
of equally spaced apertures leading into a c-on resonator volume but with 
a few partitions in the resonator cavity to prevent hot-gas flow into the 

cavity. 

With a fun length acoustic linear, the resonators are uniformly distributed 
ever the entire chamber wall from the injector to the nozzle entrance plane. 
However, shorter liners have also been proved effective. This was initially 
demonstrated by Ford (BefA-l) during a series of tests in which an effective 
(stable configuration) Ml length liner was replaced by a quarte^length 
section of itself. Stability was obtained when this section was located in 
either of the two upstream quarter-length portions of the chamber. When it 
was moved downstream to the third quarter-length position, marginal stability 
w obtained while, when located mijacent to the nozzle, hi^-frequency com- 
bustion instability was encountered. 


m addition. Garrison (Hef.A- 2 ) has reported hot firing evaluation of liners 
with lengths equal to 1/2, Vi, end 1/8 of the chamber length, but which 
were otherwise equivalent. Stability was obtained with the 1/2 and l/A 
length liners when each was installed adjacent to the injector but only mar 
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ginal stability was obtained with 1/8 length version located in this position. 
Nonetheless, the measured instability amplitudes with the l/8_length liner 
were significantly reduced from those observed with no liner (i»e^, ^ W% 
of chamber pressure peak to peak compared with ** with no liner) . 

Additional tests were also conducted with this l/i^length liner but with the 
liner moved three inches downstream. The measured amplitudes doubled with 
this change in axial position, although the amplitudes were still less than 
10 % of chamber pressure peak to peak. 

Vincent (Ref.A-3 tested 1/4, l/ 2 and. full length liners as part of a rela- 
tively extensive liner evaluation program at NASA— Lewis Research Center. 

With the liner located adjacent to the injector, stability was obtained with 
all three liner lengths. 


3n other tests at NASA— Lewis, Phillips (Ref »A-4) studied liner effects on 
stability with a liner composed of 272 individually adjustable Helmholtz 
resonators. He concluded that the damping provided by the row of resonators 
adjacent to the injector (row l) was dominant. The next two rows (rows 2 and 
3) produced negligible damping compared to the first. Hie remaining rows 
(rows 4 through 8) contributed little or no damping* 

Generally it has been concluded f*om all of the forgoing results that a 
part ia l length liner is more effective when located adjacent to the injector 
than when located downstream. The test results certainly suggest this effect; 
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however, an alt ernative ia possible: the acoustic_properties (Impedance and 

timing) of any resonator are affected by position in the engine so that the 
liner .may be less effective for that reason. The liner design should be 
adjusted for the new impedance and tuning to give comparable damping, 

Much of the foregoing work was done with uncooled liners. However, several 
programs have been conducted which were concerned largely with development 
of durable liners through the use of cooled configurations (including abla- 
tive or self cooled). Garrison (Ref. A-5) evaluated transpiration_cooled, 
f ilm cooled, and ablative liners in a series of hot firings. The transpira- 
tion cooled liner stabilized the engine when no coolant flow was used but 
became very unstable with low flow rate of hydrogen as transpiration coolant. 
Similarly, with the film cooled configuration the engine stability was 
degraded by increased coolant flow. Full stability was only obtained with 
no coolant flow. Some difficulties were encountered with the ablative liner 
as well but a stable configuration was found* 

Garrison (Ref.A-6) has also evaluated a regeneratively cooled half-length 
liner designed for use in Advanced Agena hardware. Dynamic stability was 
demonstrated with the cooled liner and also with an uncooled version of it. 
But an uncooled l/4-length liner was not sufficient to stabilize the engine. 
In earlier work (Ref. A-7) a water cooled liner was successfully used. 
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A study was recently completed at Aerojet (Ref • A— 8) during which the feasi- 
bility of acoustic liners for long duration firings was demonstrated. Two 
liners were demonstrated, one ablative and one-water-cooled configuration, 
each with total test durations of > 600 seconds. Each of these liners was 
comprised of two rows of Helmholtz resonators. The ablative liner could be 
described as a half-length liner while the water cooled liner was roughly 
one-quarter length. Hie ablative liner was found to produce dynamic sta- 
bility (but only one bomb test was made). On the other hand, the water 
cooled liner failed to damp a bomb disturbance, but adequate durability was 
demonstrated. This program included a relatively lengthy evaluation of heat 
transfer and material requirements. 

’BingHna stability has also been obtained with somewhat different absorber 
arrangements in which one or a few rows of resonators were mounted along the 
periphery of the injector. The acoustic cavity fits into this category. 

Most of the work with this absorber arrangement has been done at Rocketdyne. 
However, some success was obtained by Senneff (Ref. A-9) with Helmholtz 
resonators in the Bell Aerosystems counterparts of the Rocketdyne RS-14 and 
LM ascent engines to be described below. 

The Rocketdyne experience with acoustic absorbers of the acoustic-cavity type 
is suxonarized in Table A— 1. The absorbers used in the RS— 14, the RS— 21 and 
the 1M ascent (baffled) engines may be considered highly successful because 
the absorber completely eliminated the stability problem and no other problems 
have been encountered with it. The application to the FLQX-LPG research 


TABLE A-l HOCKETDYNE APPLICATIONS OF ACOUST 


Engine 

J RS-14 Post Boost 
Propulsion System 
Axial Engine (also 
RS-21) 

LM Ascent engine 
(with baffle) 

FLOX-LPG research 
engine 


LM Ascent engine- 
acoustic cavity 
study 

Extended Range Lex 
(XRL) Booster 

Spartan Segment 


J-2 - basic 


DTRM solid 
propellant motor 


FLEXB4 sustainer 


XRL sustainer 


Absorber 

Configuration 

Experienced 

Mode 

Design 
Mode - 

Degree of 
Success 

single axially directed 
quarter wave slots 

1 st tangential 
(9000 Hz) 

same 

excellent 

single row of quarterwave 
resonators 

3 rd tangential, 
1 st radial 

same 

excellent 

single row of slot -like 
Helmholtz resonators 

1 st tangential & 
1 st tangential- 
coupled mode 
1 st longitudinal 

same _ 

reduced peak-tc 
peak amplitudes 
from 200 $ to 30 

many quarterwave slot 
configurations (see main 
text) 

1 st tangential 
(3400 Hz) 

same 

excellent 

combination of "L" shaped 
Helmholtz resonators 

4th & 6 th 
tangential, 
1 st radial 

same 

excellent 

combination of quarterwave 
& Helmholtz resonators 

5000 & 9000 Hz 
transverse modes 

same 

incomplete - 
substantial 
improvement 

slot-type quarterwave 
resonator 

none 

1 st 

tangential 

no instability 
experienced 

85 "J"-shaped quarterwave 
tubes 

1750 Hz 1st 
tangential 

same 

failed to 

suppress 

instability 

intermediate-type resonator 
at periphery of injector 

4500 Hz injector 
coupled mode 

2 nd 

tangential 
(3300 HZ) 

failed to 
suppress - 
4500 Hz 
oscillation 

Helmholtz resonator around 
nozzle entrance 

1 st longitudinal 
(700 Hz) 

same 

significantly 

improved 

stability 

single row of quarterwave 
resonators 

1 st tangential 
(6500 Hz) 

same 

no instabili- 
ties after 
absorber 

single row of Helmholtz 
resonators 

1 st tangential 
(4500 Hz) 

same 

qualified 

success 
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enced 

e 


gentlal 

*> 


gential, 

ial 


gential & 
gential- 


gitudlnal 

gential 


rse modes 


Design 

Mode 

Degree of 
Success 

Propellant 

Combination 

Thrust 

(lb f ) 

Chamber 

Pressure 

(psia) 

Absorber Open 
Area Relative to 
Injector Area 

same 

excellent 

n 2 o 4 /mmh 

300 

125 

7$ 

same 

excellent 

NgOi^O^i*- 

50$ UDMH 

3500 (vacuum) 

122 

5 

same 

reduced peak-to- 
peak amplitudes 
from 200 $ to 30 $ 

flox(8o$ fJ/ 
55$ CH. -45$ 

C 2 H 6 

1500 to 5000 
(sea level) 

50 to 200 


same 

excellent 

H 2 0 4 / 5 0)to A - 

50^ UEMH 

3 500 (vacuum) 

122 

18 

same 

excellent 

irfna/udmh 

--- 

930 

6$ 

same 

incomplete- 

substantial 

improvement 

CTF/N g H^ 

4000 (nominal) 

500 to 700 ' 

11 & 
5.4 

1 st 

tangential 

no instability 
experienced 

k 2 V 50 ^ 2 v 

50 $ UDMH 

1600 

400 

4.2 

same 

failed to 

suppress 

instability 

L0X/GH g 

205.000 to 

230.000 

715 to 780 

4.8 

2 nd 

tangential 
(3300 Hz) 

failed to 
suppress - 
4500 Hz 
oscillation 

lox/gh 

2 

180.000 to 

265.000 

850 to 1240 

5 

same 

significantly 

improved 

stability 

ap/ctpb/ai 

* 

* 


same 

no instabili- 
ties after 
absorber 

CPF/MHF-7 

#* 

** 

8 

same 

qualified 

Success 

IRFNA-UDMH 

--- 

900 

10 


injector 


itudlnal 


gential 

) 


entlal 
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Meeting , CPIA No. 201, Vol. II, Chemical Propulsion Information Agency, Silver Spring, Maryland. 

jpjpoui FRAMt a. 








/ 


engine was not . completely successful because instability amplitudes of ~ 30 
percent of chamber pressure peak to peak were experienced with the absorber; 
however, this was a significant improvement because amplitudes ^ 200 percent 
of chamber pressure were experienced without the absorber. During the 121 
Ascent Engine Acoustic Cavity Study, which was the forerunner of the current 
program, dynamic stability was demonstrated with five somewhat different 
cavity configurations and with .an ..uribaf fladLyersion of the LM ascent engine 

injector. Eor the XRL booster engine a completely successful absorber was 

developed to_suppress several modes of instability encountered in earlier 
tests; the absorber is used in conjunction with a fouf-bladed baffle. The 
absorber used in the highly unstable Spartan sequent was partially success- _ 
ful in that the incidence rate of instability was reduced with the absorber 
but the problem was not eliminated. Hie absorber used in the RESA-5 was 
included as a precautionary measure; no instability occurred but no tests 
were made without it. The remaining engines shown in Table A— 1 were less 
than completely successful for one reason or another. Hie absorber used in 
the J-2 basic engine failed, perhaps because of inadequate open area of the 
absorber or inadequate coolant flow through the absorber; only one configura- 
tion was tested. The J2-S absorber failed to suppress an unusual "buzz" 
mode of instability with a frequency of 4500 Hz; however, that absorber was 
actually designed to suppress a bomb-induced second tangential mode (3300 Hz ) 
experienced during some tests. The absorber used in the D1RM solid-propellant 
motor gave promising results but a satisfactory design was not obtained, in 
part, because of structural failures. No instabilities have been, encountered 
with the FLEXEM engine containing an absorber but neither has the need for 
it been demonstrated. The absorber in the XRL sustainer engine has signifi- 
cantly Improved the stability of the engine but seme further improvements would 
be desirable (although not required). 
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Many of the absorbers shown in Table A-l have required no special cooling. 

The absorber designs used in the J-2 basic and J2-S engines included a 
hydrogen purge intended to control the gas properties within the absorber 
rather than cool them. The absorbers in both XRL engine are machined into 
the .ablative chamber wall material. The absorber used in the Spartan seg- 
ment required water cooling, however. 

Absorber Design Procedure'S 

Generally absorbers have been designed, from a stability standpoint, by 
selecting a configuration believed to provide sufficient damping based 
on some calculation of the damping produced by the absorber. One way of 
estimating this damping is to calculate a temporal damping coefficient, as 
described for acoustic cavities in the main text of this report. A similar 

calculation for acoustic liners is described in Reference A-10. A second 

means of estimating the damping provided by an absorber is to calculate an 
absorption coefficient, which. has been widely used for the design of liners. 
Calculation of the absorption coefficient is described by Garrison (Ref. A-ll), 
However, there has been some disagreement among various investigators con- 
cerning the usefulness of the absorption coefficient (Ref. A-12). This 
question is less important for single rows of resonators, such as acoustic 
cavities, because the absorption coefficient must be based on an arbitrary 
area in this case and, therefore, should not be used (c.f. Ref. A-13). 
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Priem and Rice (Ref. A-U) and Phillips (Ref. A-15) have analyzed the damping 
produced by an acoustic liner. Both analyses resemble, to some extent at 
least, the temporal damping coefficient calculations . Zinn (Ref. A-16), 
Sirignano (Ref. A-17) and Smith, et. al. (Ref. A-18) suggest that liners be 
designed by maximizing the real part of the liner acoustic admittance. How- 
ever, this recommendation results from considerations similar to, but more 
restrictive than, the damping coefficient calculation. 


ABSORBER IMPEDANCE 

The acoustic impedance of the absorber elements may be written as 

z a * R + jX * A_1 

where R is the resistance and X 13 the reactance. Each of these components 
has been the subject to considerable study; however, a complete understanding 
of their behavior under motor conditions is not yet available. Each is 
affected by the steady, flow conditions, the temperature distribution and the 
local pressure amplitude. However, the amplitude dependence of the reactance 
term is small and it is adequately described, at least in form, by a linear 
(lew-amplitude) expression. 

The impedance for a Helmholtz resonator may be well approximated as 

» z. ♦ z„ A -2 

a A c 

where z is an impedance for the aperture and z Q is an impedance for the 
£\ 

resonator cavity. 
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Cavity Impedance 


The Impedance of a simple resonator cavity is well known and is given by 



1 P 0 s 

wV 


A- 3 


where V is the volume o± the cavity, u> is the angular frequency, S is the 
cross-sectional area of the aperture, and c 2 = 7 p/p is the isentropic 
sound velocity in the cavity. Equation A - 3 applies if all dimensions of the 
cavity are small compared to a wavelength (Based on the local sound velocity). 
If. the depth of the cavity does not meet this restriction, and the cavity 
is cylindrical or rectangular in shape, Eq. A - 3 may be replaced by 

z = -j JL2§ cot <t2k A- 4 

c A — c 

where A is the cross-sectional area of the cavity and L is the cavity depth. 
Equation A- 4 is readily obtained by solving the wave equation for the cavity. 
For small arguments, Eq. A- 4 is equivalent to Eq. A- 3 . Equation A- 4 was 
obtained by vising a short wavelength approximation. 

Aperture Reactance 

The reactance of the aperture in a Helmholtz resonator which is short com- 
pared to a wave length may be written as 

X A = P A w t Q A - 5 

where a * + 8 , the effective aperture length j P A is the time-averaged 
density in the aperture} and u is the angular frequency of the oscillation. 
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The effective aperture length is obtained by adding a mass end correction, 5 , 
to the physical length. 

The mass end correction depends on both amplitude and steady flow velocities j 
for low amplitudes and with no flow, this end correction is usually written as 

6 L b °- 85 d a (1 - 0.7/aj) A-6 

l 

Equation A-6. is based on the analysis of Ingard (Ref. A-l9)and should be 
restricted to 0^ < 0.16. 

At sufficiently high amplitudes and/or with high through- flow velocities the 
end correction diminishes to a value near half the linear value (Ref. A- 20 
through A-22). Based on the theoretical work of Westervelt (Ref. A-23), a 
limiting value of 3/8 6l is often referenced. Some confusion has pre- 
viously arisen in application of the Westervelt factor, because of the failure 
to distinguish between sharp-edged and finite-length orifices. That is, at 
times the factor 3/8 has been applied to the linear effective length for a 
finite length orifice rather than only to the linear end correction, this 
correction being appropriate for sharp-edges orifices only. 

The effect of steady cross flow is less clear; however, at least qualitatively, 
it appears the cross flow tends to sweep away the normal end correction on 
the side with flow*- . If the outside end correction is completely swept away 
by the steady flow, the end correction should be 

6 cp * 0.8S D a (0.5 - 0J7y/al) A-7 
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(the correction for pertains only to the cavity-side end correction). 
Phillips (Ref . A-22) suggests the use of 


s cf * i \ 


A- 8 


which is. reasonably compatible with Eq, A-7* Further, Phillips has found 
the commonly used correction (Ref.A-22), which is based on the work of Mechel, 
et al (Ref.A-2), to be excessive. 


An end correction, for high amplitudes with steady through and/or cross 
flow, of 


6 = 0.8S D a (0.5 - 0.7/?') 


A-9 


■4 

% 


appears to be a reasonable representation of the latest results. 

For noncircular apertures, a similar expression should apply* It appears 
Eq. A-9 should be generalized to 


6 = 0.96 (S A ) 1/2 (0.5 - 0.7 /op 0.16 


A- 10 


M 


I 
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Absorber Reactance 


The reactance of an ordinary Helmholtz resonator is (the sum of the aperture 
and cavity reactances) 


where 


Y p S 

x = p A w h " 



A-ll 


The reactance of a quarterwave resonator, for a partitioned cavity, is 

ui£ 

X * -p A c A cot A-12 

A 


where i is the effective length of the slot. 


Aperture Resistance 

The resistance of absorber elements is substantially affected by the local 
pressure amplitude, and by the steady-flow environment. In the amplitude 
range of Interest for absorber design, the resistance of Helmholtz and quar- 
terwave resonators is dominated by high amplitude effects, the viscous losses 
being generally negligible in comparison. 


Much of the early work on acoustic liners was based on the work of Ingard. 
Ingard measured the aperture impedance of an orifice in the nonlinear range. 
He chose to correlate his data in terns of nonlinear end corrections and 
found the resistance end correction,, A^/D , varied with velocity amplitude 
in the aperture or the displacement amplitude, depending upon whether or not 
the displacement amplitude was comparable to the aperture length. Similar 
effects were found with the nonlinear' mass end correction. Ingard^ data 
were reasonably well correlated in this fashion. He attributed the observed 
nonlinear effects to turbulence. Ke did not, however, develop an analytical 
model fOrJthese effects, and his correlation was simply a means of repre- 
senting the data. He did not show the correlation to be valid for conditions 
otherJbhan-thoseLJSxisting in his apparatus (air at ambient conditions). 

Later work by Blackman (Ref.A-24)on Helmholtz resonators and by Garrison 
(Ref. A-2,A-25,A-26) has been based on the Ingard work. Furthermore, most 
existing liner designs have been developed on this basis. However, more 
recent work clearly shows a more appropriate formulation. 

Recently Ingard has re-examined the nonlinear problem (Ref.A-20). He has 
carried out additional experiments using a somewhat different technique. He 
was able to show good correlation between his latest data, his earlier data, 
the results of Sivian (Ref,A-27)and those of Bies and Wilson (Ref. A-20). 
Furthermore, all of these data for the aperture resistance, R , are well 
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described at high amplitudes by the expression 

R = P A u A -13 

where u is the velocity amplitude in the aperture. An equation of this 
form was originally proposed by Sivian (Ref • A-27 ) in 1935 * Ingard (as did 
Sivian) indicates this result, Eq. A-13, might be expected on the basis of 
Bernoulli's, equation. Phillips (Ref.A- 29 ) has experimentally confirmed Eq. A-13. 

2 inn (Ref .A-16) and Sirignano (Ref.A-17) have examined the problem analytically. 
Their results agree with each other except for a numerical coefficient, and 
they agree with the form of Eq. A-13. Hie same form can be obtained by a 
quasisteady approximation and the equation for steady flow through an ori- 
fice (Ref. A- 35) Tonon (Ref. A-36)has recently reported an analytical treat- 
ment of Helmholtz resonators. 

It is lik ely that a discharge coefficient for oscillatory flow should be 
included in equations of the fom of Eq, A- 13 * Garrison (Ref. A-ll) has ob- 
tained a good correlation of this acoustic resistance data by introducing 
the steady flow value in the fom 

R » 0.37 ptyc£ A-Il+ 

Others have suggested the use of a discharge coefficient as well. 
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The high amplitude resistance of quart erwave slots has been measured at 
Rocketdyne (Ref. A- 30 ). The resistance results were correlated with an ex- 
pr ess ion si milar to Eq. A- 13 * i.©. ... 

R = r p A fi A-15 

where, for the slot, r » 0.6 was obtained. Thus, Eq. A-15 appears to be 
the best available expression for the nonlinear resistance of absorbers, 
where r = 0 . 37 /C D 2 for HeMioltz resonators. 

Steady_ through flow has a similar effect on resistance. Ingard (Ref. A- 20) 
also measured the orifice impedance for this condition. His results are 
given, at high velocities, by the expression 

R = P A U A-l6 


where U is the steady flow velocity in the aperture, and u << U . 

3 

Garrison (Ref. A-ll) suggests R ” p.u/C based on his data, which differs 

JrL JJ 

somewhat from, his earlier expression of Ref. A- 3 4 . 

Ingard (Ref . A-21) suggests the following approximate equation for use when 
u ~ U. 


R “ P A 0 4 »A U 


A-17 


This equation limits properly Phillips (Ref.A-22) suggests an expression of 
the same form but with the steady-flow term multiplied by 1.5* 
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The effects of steady cross flow (perpendicular to the axis of the aperture) 
are even less clear. Most liner designs have been based on the work of 
Meyer, Mechel, et al (Ref. A~31and JW32). They found both components of the 
impedance to be affected by the flow; the aperture resistance was found to 
increase with cross-flow velocity. An interesting feature of this work is 
that a m plif ication was encountered In some regions of operation, but this 
does not appear to be a problem in absorber work (Ref. A- 33). 


Recent work by Garrison (Ref .A- 11) led to the following empirical expression: 

R = R 0 (1 + 1.9 M p ) A-18 

where R 0 is the no-flow resistance. This differs from an earlier expression 
reported by him (Ref.A-25). Phillip's results (Ref.A-29) appear to obey a 
similar relationship with a somewhat higher slope. 

Garrison (Ref. A-34) also measured the aperture resistance with simultaneous 
through and cross flows. He concluded the resistance was dominated by the 

through-flow effect and was adequately described by the through-flow ex- 

pression, except for very low through-flow velocities. 


The effects of steady flow on the resistance of quart erwave resonators have not 
been. measured. However, it appears reasonable to assume the effects are 
similar to those found with Helmholtz resonators. 


Hot-Firing Impedance Measurements 


During several programs at Rocketdyne (Ref. A- 10 and. A- 30 ) and at NASA- 
Lewis (Ref. A - 4 and-A= 22 ), and recently at Pratt- and Whitney (Ref.A-13), 
oscillatory pressure measurements were made In the acoustic absorber to 
determine, insofar as possible, the acoustic impedance that existed during 
test firings of an engine. The technique is essentially that described by 
Sivian (Ref. A- 27). Related measurements were made by Crocco, et al (Ref. A- 37) 
Results from these measurements tend to be widely scattered because of diffi- 
culties in making the appropriate measurements, although Garrison's results 
show little .scatter (Ref.A-13). However, the technique is regarded as highly 

worthwhile because it provides information on the actual behavior of the - 

absorber during use; otherwise, cold-flow impedance data must be extrapolated. 

With this technique, oscillatory pressure measurements are made et the open 
end and at the closed end of the absorber. These data are used to infer the 
absorber impedance. 

Absorber Temperature and Gas Sampling Measurements 

The acoustic impedance, and, in turn, the damping of an absorber is often 
strongly dependent on the time-averaged temperature distribution within the 
absorber. Considerable temperature data are now available; from these, some 
conclusions can be drawn. 
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The variations in measured temperatures are often rather large, indicating — 
these temperatures are far from uniform and that they are difficult to 
measure. With acoustic liners, the temperature varies with open area frac- 
tion. With ouarterw ave res on ators , often, but not always, there is a strong 
temperature gradient from the open to closed end of the cavity. And, as 
might be expected, the temperature varies with the propellant combination. 


Temperature data from hydrogen/oxygen acoustic liner test firings may be 
found in Ref.A-4,A-33, and A-24. Similar data from LCX/RP-1 firings may be 
found in Ref .A- 38 andA-39. Similar data from NgO^/N^-UDWH (50-50) firings 
may be found in Ref. A- 2 ,A-22, and A-30. 


Temperature data from quarterwave-type absorbers with the NgO^/N^-UDMH 
(50-50) propellant combination may be found in Ref. A-30. 

g ome degree of temperature control can be exercised by purging the absorber 
with cold gas (Ref. A-22. 

Gas sampling has been used to determine the gas composition within the 
absorber, which is important for design calculations. Garrison (Ref. a- 5). 
used a sampling technique for ablative liners. Sampling measurements were 
made in the Rocketdyne F-l and LM Aseent absorbers (Ref. A- 38 and A-40). 
Phillips (Ref .A- 4) employed gas sampling in his work with the Hg/C^ pro- 
pellant combination. In each of these_eases, the gas composition was found 
to be significantly different from that of the equilibrium combustion 
products. 
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Pressure probes were used in the F-l engine to measure the local gas 
velocity (Ref. A-38). 

Partitioned Absorbers 

Partitions are commonly used in the resonator cavity of acoustic liners to 
support the liner and to prevent hot-gas circulation into the cavity. 

Phillips (Ref.A-22) found the addition of three circumferentially directed 
partitions reduced the measured cavity temperature. 

Partitions are not needed for acoustic or damping reasons, if appropriate 
impedance expressions are used. Partitions are desirable in slot- type ab- 
sorbers to prevent overheating and erosion. 

Design Amplitude 

Some uncertainty and confusion exists relative to the amplitude for which an 
absorber should be designed. Because of the impedance nonlinearities, the 
predicted damping varies with amplitude. Further, should the absorber be 
designed to damp (be most effective) a fully developed instability or, on 
the other hand, damp the oscillation before it develops? 

Most early liner calculations to date have been made with an arbitrarily 
assumed ’’incident" amplitude of 190 db. The "incident" amplitude was assumed 
to be half the local amplitude; thus, the 190-db level corresponds to a 
local amplitude of 56 psi peak-to-peak. 
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Another amplitude which is believed to be more appropriate, is a peak-to-peak_ . 
oscillation of 20 percent of chamber pressure (Ref. A-10 and A-30). This 
amplitude is based on the observation that the acoustic^equations in use for 
all damping calculations are not valid for higher amplitudes (Ref. A-4l). 

Further, it is the fractional amplitude rather than the absolute amplitude 
that appears in the damping calculations. 

In most cases an engine is regarded as stable if the peak-to-peak amplitudes . 
are less than lCLpercent and, further, the combustion noise levels are at 
least a few percent. Thus, it appears inappropriate to design an absorber 
for peak-to-peak amplitudes less than a few percent or greater than 20 per- 
cent of chamber pressure. 

If 190-db level is arbitrarily employed for all chamber pressures, at low 
chamber pressure the equations are invalid whereas at high chamber pressure 
the amplitude is below the noise level of the engine. 

Bandwidth 

In the usual sense, bandwidth refers to the frequency range over which an absorber 
is effective. A related characteristic is the temperature range over which 
an absorber is effective. Generally, it is considered desirable that an 
absorber design have sufficient bandwidth to damp all instability modes that 
might occur. Similarly, it is desirable the absorber have sufficient damping 
over the entire temperature range_the absorber might experience. However, 
there is some danger that overall damping will be sacrificed in favor of band- 
width to the extent the absorber will not be adequate in any range. 


The most effective means of obtaining increased bandwidth is to employ 
resonators with several resonant frequencies (Ref. A-34). With Helmholtz 
resonators, the bandwidth (and temperature sensitivity) may be improved by 
decreasing the aperture length. 

Cooling and Materials Selection 

There are few special cooling and materials selection requirements that are 
unique to acoustic absorbers. Generally, standard thrust chamber design 
techniques are effectively applied. 

Often, no special cooling provisions are required for single -row absorbers 
mounted in the injector. In one case, a thin layer of ablative is used on 
the combustion side* 

Studies of cooling techniques for acoustic liners are described in Ref.A-2 , 
A-7 ,A-25,A-34,A-39, and A-42. Garrison describes the principal experience 
with cooling and thermal protection of liners in Ref.A-4-3, Several water- 
cooled liners have been successfully tested. 

There are some stability-related problems associated with film eooled and 
transpiration cooled liners. According to Ref.A-25,the coolant flow tended 
to decrease the stability of the engine. This effect is not understood but 
may be associated with a convection of acouster energy into the system due 
to the coolant flow (Ref. A-4l). 
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NOMENCLATURE FOR 



APPENDIX A 

absorption coefficient 
isentroplc sound velocity in chamber 
isentropic sound velocity in aperture 
discharge coefficient 
aperture diameter 

(-d" ! 

cavity or aperture length (depth) 
effective aperture length 
resonator cavity length 
Mach number of steady cross flow 
time average pressure 
resistance 

aperture cross-sectional area 
amplitude of oscillatory velocity 
steady flow velocity 
resonator volume 
reactance 

aperture impedance 
resonator cavity impedance 
heat capacity ratio 
resistance-coefficient 
mass end correction 
with cross flow 

linear (low amplitude) mass end correction 
steady gas density in chamber 
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steady gas density in aperture 
fractional open area from resonator side 
angular frequency 
resonant angular frequency 
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APPENDIX B CAVITY TEMPERATURE DISTRIBUTIONS * 

Calculation of the damping characteristics of acoustic slots requires_that 
the gas temperature distribution in the cavities be described with reasonable 
accuracy. The analysis presented in this appendix was made to develop a con- 
sistent method of predicting temperature distribution in acoustic cavities as 
a function of cavity dimensions, cavity location, and propellant combination. 

As might be expected, the results of the analysis suggest that the injector 
configuration is also a very important parameter, the effects of which cannot 
as yet be described in a general manner. 

The temperature distribution within a cavity must be a function of the heat 
transfer characteristics of the flow field within the cavity. The analysis 
was begun by reviewing the literature concerning flow and heat transfer within 
cavities and slots; the findings of this review are discussed in the following 
paragraphs. In general, the data available in the literature were of limited 
direct application to the present problem because they involved well-defined 
external gas flow fields rather than the accelerating, three-dimensional, 
two-phase flow field which prevails adjacent to a liquid rocket injector. The 
literature survey is followed by derivation of a simplified model for the 
cavity t emp erature distribution which is then correlated with the experimental 
data obtained during the full scale hot firings. 


*Work done by W. S. Hines of Roeketdyne. 
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LITERATURE SURVEY 


The most pertinent reported investigations of flow and heat transfer in rec- 
tangular cavities were those of Fox (Ref. 3-1), Haugen and Dhanek (Ref. B-2, 
B-3) and Burggrof (Ref. B-4). Other_jcaported studies ordinarily dealt with 
high velocity flows and involved notches with aspect ratios (b/H) greater 
than 1.5; therefore, their results were expected to have no more than limited 
applicability to the present problem which involves aspect ratios less than 
unity. 

Fox (Ref. 2) studied the heat transfer within notches of b/H between 0.25 and 
1.75 with uniform free stream velocities ranging from 160 to 580 ft/sec. The 
local heat transfer coefficient along the cavity walls (z)* defined as 



was found to very as u 0 and z”° along the downstream wall of the cavity. 
The cavity gas-temperature distributions were not measured at low aspect ratios 
because of experimental difficulties. However, the cavity gas* temperature 
results for b/H « 1.0 Indicate that the cavity, gas temperature is closer to 
the main stream gas temperature than to the_wall temperature (even at the 
bottom of the cavity) which implies that most of the resistance to heat trans- 
fer is in a film at the wall. 


* 

Coordinate z is measured from the down stream lip of the cavity with the 
positive z direction taken into the cavity. 
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Haugen (Ref. B-2, B-3) analyzed turbulent flow and heat transfer in rectangular 
cavities. He supported his analy sis with measurements of pressure recovery 
coefficients, of heat transfer coefficients along the cavity wall and of gas 
t emp era ture pro files within the cavity. Haugen also performed flow visualiza- 
tion experiments which indicated that the flow within the cavity has a cellular 
eddy structure of the type shown in Fig, B-l. Because of the presence of 
these eddies., the temperature distribution has a stepwise variation with cavity 
depth as shown in Fig. 2 , which was taken from Ref. B— 2. Haugen only reported 
his calculated velocity and temperature distributions to depths (y/b) less than 
1.0 j however, to those depths, the agreement between experimental and calculated 
distributions was excellent. Haugen showed both analytically and experimentally 
that the flow and heat transfer within the cavity were functions of the boundary 
layer thickness upstream of the cavity, i.e. of the main stream flow. 

Burggrof (Ref. B-4) analyzed the velocity and temperature distributions in a 
square cross section cavity. A modification of his egression, 

T - T , A = 2 laminar . . 

_ 1 (B-2) 

T ® " T w ~ n/A( 1 + b/H . _ + 

® w A = 3 turbulent 


was used by Haugen to define the average temperature in the "in viscid core'* of 
the cavity defined by the counter-rotating eddies of Fig. B— 1. 

The cavity gas temperature defined by Eq. B=2 is nearly uniform. This does 
not agree with the experimental temperature distribution obtained by Haugen as 
shown in Fig, B-2. It definitely dose not agree with the acoustic cavity 
temperature distribution data from the present program as exemplified in Fig. B-3. 
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Figure B-2 TEMPERATURE DISTRIBUTION WITHIN A RECTANGULAR 
SLOT FROM HAUGEN'S EXPERIMENTS (Ref. B-3) 
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The probable reasons for the lack of agreement with the acoustic cavity data 
are the previously mentioned differences in the main f l ow field conditions* 
However* even for a well defined external flow fi eld* Eq. B-2 cann ot be a pplied 
for very deep cavities , since it predicts a limiting value of A ^ for the 
dimensionless temperature rather than unity. 

ACOUSTIC CAVXTI TEMPERATURE DISTRIBUTION 

The purpose of. this section is to describe the derivation of a relation for 
temperature within an acoustic cavity which provides a simple, yet plausible, 
correlation for the experimental temperature profiles shown in Fig. B-3 and 
which can be adapted for correlation of more extensive temperature data when 
available. 

Because the acoustic slot is located next to the injector face where the 
degree of combustion is unknown, there is actually little that is well defined 
about the problem except for the cavity dimensions. Therefore, a simplified 
model was used wherein only one— dimensional heat conduction was considered but 
the coefficients were treated as empirical parameters. For the cavity con- 
figuration shown in Fig. B-4, the thexmal energy equation is 

dT 

|y ( 'Vc “ay = ph( V T v } + Across (B-3) 


t 

I 
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Figure B-3 DEFINITION OF HEAT TRANSFER SYSTEM WITHIN AN ACOUSTIC SLOT 
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Figure B-4 MODEL FOR VOLUMETRIC HEAT GENERATION TERM IN ACOUSTIC 
CAVITY BASED UPON EVAPORATION AND DECOMPOSITION OF 
ENTRAINED FUEL FILM COOLANT 


a 
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n 




where 


k T (y) is the total thermal conductivity for heat transfer through 
the cavi ty g as which is assumed to include molecular thermal con- 
ductivity, the eddy conductivity and any recirculating cavity flow. 


Q(y) is the volumetric heat generation in the cavity gas which is 
assumed to include combustion, monopropellant decomposition, and 
spray evaporation. It may, therefore, be positive or negative. 


Employing two simplifying assumption that 


d k-Xy) 


km(y) dy 


* constant * - a or k_ = L e 


•ay 


cross T 


h g (y) 

sfer 


constant ■ (3 


which do not appear unreasonable (at least for limited ranges of y) from the 
usual convection theory. 

The differential equation may be written as 


- P9 - <Hy) 


dy 


,2 dy 
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where 

0 " T g - T W (H) 
T y (y) * AT y (y) + T y (H) 

°(y) “ - fwr(y) 

and T (H) is assumed constant. 

V 


The coefficients a and (3 and the "catch-all" function G are yet to be 
determined. Solution of the homogeneous equation is 


where 


-c 2 y 

0 8 o^e + c^e 


e = - a _ Va 2 + 1(3 
2 2 2 



(B-5) 


Because 6 must be bounded, 
then 


» 0. 


Thus, with no heat source (G ** 0), 


6 = 




(6-6) 


Haugen (Hef. 4) separa ted the flow in and around a cavity into three zones: 
(a) the free stream outside the cavity (b) the mixing zone or shear layer 
and (c) the inviscid cavity flow. The eddy diffusivity e was different in 
each zone. Although a well defined shear layer probably does not exist 


immediately downstream of an Injector, the acoustic cavity temperature data 
suggest that at least two zones should be considered for the cavity flow 
field: 

For zone 1 (o < y < y x ) 

-c 0 y 

e * e 0 e 

For zone 2 , . 

(y x < y) 

-«a(y-yi) -( c 0 “ c 2^ y i " c 2y 

e ® e^e * e Q e e 

Haugen defined the thickness of the shear layer in terms of a free pla. ar jet 

y l • * o88x t -9) 

In the present problem, where a well defined free stream velocity does not 
exist and the temperature field is considered one-dimensional, the boundary 
between zones 1 and 2 will be defined, for convenience, as 

y i = b/2 

This corresponds approximately to one half the size of the c ellul a r eddies of 
Fig, B-l, i.e., to the maximum plausible size of any layer corresponding to a 
shear layer. 


r: 


Examination of Eq. B-7, B-8, and B-9 together with the data shown in Fig. B-3 
demonstrate the need for heat generation term, Q. Any relation based strictly 
upon one-dime nsional transport of sensible energy from the main flow will 
result In prediction of a monotonic temperature decrease through the cavity, 
a pr ediction which is at v ariance with the data. 


The physical model to be employed in the. calculation for the heat generation 
term Q assumes that a small amount of fuel film codai.. or propellant 
spray is deflected into the cavity where it evaporates -and then reacts either 
indothermically or e xothe rmically (Fig. B-5). Because the dependence of Q 
on y is unknown except that it must approach zero as y becomes large , the 
following simple function is arbitrarily selected: 


o-ve-v. 


The corresponding solution of the differential equation 


-c 0 y -CKy 

6 - e e + o-e y < b/2 


(B-10-a) 


e= e 0 e 


-( c o - c 2 )b/2 -c 2 y 


: - c 2* -<V 

e + e e y > b/2 


(B-10-b) 


where the coefficients and differ from those used previously (and 
then eliminated). 
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The form of Eq. B-10 is well suited for simple correlation with limited 
experimental cavity temperature data, although it cannot provide as accurate 
(or theoretically sound) a correlation as more complicated heat generation 
expressions, e.g., 

•\y s -*i 2 y 

Q s Y i e 1 + V g e 


Q = v 3 y 


e"V 


which functionally permit the heat generation to exhibit a maximum. Similarly, 
when_sufficient data are available, the analysis should be made two-dimensional 
within the cavity. 


CORRELATION OF EXPERIMENTAL DATA 

The constants in the cavity temperature distribution relation (Eq. 8a and 8b) 
were ob tain ed by curve fitting the temperature data obtained from the short 
duration motor firings described earlier in this report. A variety of cavity 
configurations were tested. Temperatures were measured with tungsten rhenium 
thermocouple probes placed in the cavity. Up to six local measurements of 
temperature were made within a cavity at depths of 0,1-, 0.2-, 0.3 -, 0,6-, 
0.9-, and 1.2-inch from the aperture opening in the chamber wall* Measured 
temperatures for the radial cavity configurations of Table 7 are listed in 
Table B-l, 


214 


ro ro ro 


TABLE B-l 














A number of factors limited the effectiveness of the correlation between 
Eq. B-10 and the data in Table B-l. First, steady state was_aasumed while the 
data were obtained firm short duration firings and, therefore, must represent 
a transient heat transfer situation. Because the probe temperatures leveled 
off at essentially constant values in a fraction of a second, the assumption 
was made that the heat transfer coefficient between the cavity gas and the 
cavity walls was much lower than the coefficient between the cavity gas and 
the main stream and, moreover, was low in absolute magnitude. This implies 
a quasisteady condition for the cavity gas and also implies a cavity wall 
temperature near the original value} a cavity wall temperature of 100 °F was, 
therefore, assumed for use in the parameter. Secondly, the thermocouple 
probes were all inserted the same distance from the upstream cavity walls . 

Any two— d i -i nwa.1 effects associated with the type of eddy structures shown 
in Fig. B-l are, therefore, not fully described. Finally, the temperature of 
the combustion gas and the spray concentration in the main chamber immediately 
outside the aperture opening are also unknown. 

Because of these factors, a number of simplifications were used in the data 
correlation. The gas temperature at y = 0 (the aperture opening) was ob- 
served to be 5200 o F for all the cavities, i.e., at approximately the thermo- 
dynamic flame temperature for the bulk: mixture ratio. This specifies e Q of 
Eq. B-10 to be 5100°F. The heat generation turn was assumed to be essentially 
aero for y greater than b; therefore, the coefficient C 2 was evaluated from 
the slope of the best-fit straight line on semilogorithaic coordinates through 
the experimental e(y) for y greater than b. Because the range of the ex- 
perimental data was. limited, the correlation was based upon an arbitrary single 

*tlpstream is definod as being the side of the cavity closest to the injector. 
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value of An Indication of the wall temperature was obtained from the 
erosion characteristics of an ablative chamber using the same injector. These 
characteristics indicate that the gas temperature at the wall was at least 
equal to the melting temperature of the phenolic refrasil (approximately 
3400 °F)j therefore, was taken to be 

e 0 + <= 3 = 33oo°f 

= -i8oo°f 

lhe coefficients and were then chosen from the temperature plots 
(semi logarithmic) to give a satisfactory fit to the experimental data. 

Hie resultant coefficients are listed in Table B-2. The corresponding tempera- 
ture distributions are compared with experimental temperatures in Fig. B-6 - 
B-13. As might be expected, c 2 and c^ vary inversely with cavity length b, 

because they should be approximately inversely proportional to Across which 

P 

is nearly proportional to b. Similarly, the decay of the temperatures with 
increasing depth (larger values of and ) becomes greater with in- 
creasing distance from the injector* Hiis is expected because the main com- 
bustion gas flow outside the cavity more closely resembles one-dimensional 
stream with a thinner shear layer at the cavity opening. 

Examination of Fig* B-6 - B-I3 shows that Eq. B-10 represents the data fairly 
well at y > b , but the agreement is poor for y < b. Hie agreement is par- 

ticularly poor for the thermocouple at y = 0,1-inch. Better agreement could 
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Figure B-5 CORRELATION OF TEMPERATURE DISTRIBUTION IN ACOUSTIC 
CAVITY NUMBER 4 WITH EQUATION B-10 
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Figure B-6 CORRELATION OF TEMPERATURE DISTRIBUTION IN ACOUSTIC 
CAVITY NUMBER 5 WITH EQUATION B-10 
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Figure B-9 CORRELATION OF ACOUSTIC CAVITY • TEMPERATURE DISTRIBUTION 
IN ACOUSTIC CAVITY NUMBER 14 WOT EQUATION B-10 
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FlKure B-10 CORRELATION OF TEMPERATURE DlfflRIBUTION IN ACOUSTIC 
CAVITY NUMBER 9 WITH EQUATION 1-10 
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Figure B-ll CORRELATION OF TEMPERATURE DISTRIBUTION . IN ACOUSTIC 
CAVITY NUMBER 10 WITH EQUATION B-10 
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Figure B-12 CORRELATION OF TEMPERATURE DISTRIBUTION IN ACOUSTIC 
CAVITY NUMBER II WITH EQUATION B-10 


I 


be obtained by manipulating coefficients c Q and c 4 in Eq. B-10-9; however, 
it is likely that the temperature indicated by the thermocouple at y “0.1- 
inch is lower than the desired mean temperature across the cavity cross 
section at this depth. The thermocouples were introduced into the cavity from 
the upstream wall with the exposed junctions uniformly located a distance of 
0. 2-inch from the upstream wall. Because of the relatively high conductivity 
of the thermocouple material, the junction should indicate a mean temperature 
for the cross section from the tip back to the upstream wall. This mean 
temperature will be substantially lower than the desired average temperature, 
because of the circulation patterns in the cavity, for the overall cavity 
cross section at the given value of y. According to the temperature data 
obtained by Fox (Ref. B-2), this variance may be particularly significant in 
the region just beneath the upstream life of the cavity which is near the 
thermocouple at y = 0.1*. The effect of distance from the upstream wall is 
due to the mixing layer, which has a thickness proportional to the distance 
from the upstream wall, as defined by Eq* B-9* 


A reasonable method is needed for predicting the approximate values of the 
coefficients 0 o , C 2 , and in terms of the flow and geometry parameters 
which are expected to define them. FronuEq. B-5 the coefficients C 0 and C 2 
are expected to depend on the cross-sectional area of the slot. In fact, 
when k T is independent of cavity depth y , C q and C 2 should be inversely 
proportional to the slot width b. Similarly, because the thickness of the 


*In Fox’s measurements, temperatures were high near the upstream wall 
because the experiments represent heat transfer from a heated wall to a cold 
gas. The direction of the heat transfer is the opposite in the acoustic 

cavity. 
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mixing or shear layer at the cavity inlet should be decreased by the presence 

of a well-defined axial flow and chamber wall boundary layer, the coefficients 

should increase with distance from the injector. Fig. B-13&, b, c show the 

variation of the coefficients C q , C^, and with slot width while Fig. B-l6a, 

b, c shows their variation with distance from the injector. The parametric 

variations are approximately as expected. The variation in coefficient C, 

4 

with slot depth and location is similar to that for C and <L. 

o 2 

The curves shown in Fig. B-13 and B-14 are not, by themselves, completely 
sufficient to permit accurate prediction of cavity temperature distribution 
because they may be specific to the given injector and they do not satisfy 
the above-mentioned problem of the possible two-dimensional temperature dis- 
tribution in a given cavity. However, these correlations provide a guide- 
line for estimating the temperature distributions. 
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nomenclature for appendix b 


X 

A 


b 

C o ,C l ,C 2 ,C 3 ,C ^ 

G 


h 

g 

H 


M 

P 

q. 


coefficient in Eq.. 2 
area 

cross sectional area of cavity 

cavity length 

coefficients 

general function 

heat transfer coefficient 

cavity depth 

turbulent conductivity 

coefficient 

perimeter 

heat flow 




heat generation term 


temperature 

distance in direction of gas flow 
distance into cavity 
distance along wall from downstream 
lumped coefficients in Eq.. B-4 


Subscripts 


gas 

wall 

at y * 0 
mainstream 
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